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Scilab numbering policy used in this document and the relation to the 
above book. 

Exa Example (Solved example) 

Eqn Equation (Particular equation of the above book) 

AP Appendix to Example (Scilab Code that is an Appednix to a particular 
Example of the above book) 

For example, Exa 3.51 means solved example 3.51 of this book. Sec 2.3 means 
a scilab code whose theory is explained in Section 2.3 of the book. 
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Chapter 2 

Steady State conduction In one 
dimension 


Scilab code Exa 2.1 STEADY STATE RATE OF HEAT GAIN 

//Example 2.1 

//(a)calculate the steady state rate of heat gain 
//(b), th etemp. of interfaces of composite wall. 
//(c)the percentage of total heat transfer 
resistance . 


/ / additional 
// Given 

thickness of cork. 


A = 1 

//m" 2 , area 


//for inner 

layer (cement) 


ti=0.06 

/ /m, thickness 


CM 

o 

ii 

■H 

/ /W/m C, t h e r m a 1 

conductivity 

Ti=-15 

//C, temprature 


//for middle 

layer (cork) 


tm=0.1 

// m, thickness 


km=0.043 

//W/m C, thermal 

conductivity 

//for outer 

layer (brick ) 


to =0.25 

/ / m, thickness 


ko =0.7 

/ /W/m C, thermal 

conductivity 

o 

CO 

ii 
o 
E— 1 

/ / C, temprature 
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/ / Calculation 

//Thermal resistance of outer layer //C/W 

Ro=to/(ko*A) 

//Thermal resistance of middle layer //C/W 

Rm=tm/(km*A) 

//Thermal resistance of inner layer //C/W 

Ri=ti/(ki*A) 

Rt=Ro+Rm+Ri 

tdf=To-Ti //temp driving force 

//(a) 

Q=tdf/Rt //rate of heat gain 

printf(”the rate of heat gain is %fW\n”,Q) 


//(b) 

// from fig. 2.4 
tdl=Q*to/(ko*A) 
layer 
T1=To-tdl 
and cork 
II similarly 
td2=Q*tm/(km*A) 
layer 
T2=Tl-td2 

cement and cork 
printf (”interface 
C\n”,T1) 

printf (”interface 
%f C\n”,T2) 


//C temp, drop across the brick 
//interface temp, between brick 

//C temp, drop across the cork 
//C, interface temp, between 
temp, between brick and cork is 
temp, between cement and cork is 


//(c) 


Rpo=Ro/Rt 

II thermal 

resistance 

offered 

by 

brick layer 

Rpm=Rm/Rt 

II thermal 

resistance 

offered 

by 

cork layer 

Rpi=Ri/Rt 

II thermal 

resistance 

offered 

by 


cement layer 
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printf (” thermal resistance offered by brick layer is 
%f percent\n” ,Rpo* 100) 

printf (” thermal resistance offered by cork layer is 
%f percent\n” ,Rpm* 100) 

printf (” thermal resistance offered by cement layer 
is %f percent\n” ,Rpi* 100) 


// second part 

x = 30 

qi=q*(i-x/ioo) 

Rth=tdf/q1 
Rad=Rth-Rt 
Tad=Rad*km*A 

printf (” Additional thickness of cork to be provided 
=%f cm”,Tad*100) 


//percentage dec in heat transfer 
//W, desired rate of heat flow 
//C/W, required thermal resistance 
/ / additional thermal resistance 


Scilab code Exa 2.2 Rate of heat loss 


//Exm [pie 2.2 
//Page no. 15 
// Given 

//outer thickness of brickwork (to) & inner 
thickness (ti) 

t o = 0.15 // m 

ti=0.012 //m 

//thickness of intermediate layer(til) 

til=0.07 //m 

//thermal conductivities of brick and wood 

kb=0.70 //W/m celcius 

kw=0.18 //W/m celcius 

//temp, of outside and inside wall 

To=-15 //celcius 

Ti=21 //celcius 

//area 

A = 1 / /nr" 2 
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//(a)solution 

//Thermal resistance of brick , wood and insulating 


layer 

TRb=to/(kb*A) //C/W 
TRw=ti/(kw*A) //C/W 
TRi=2* TRb //C/W 

//Total thermal resistance 
TR=TRb+TRw+TRi //C/W 
//Temp, driving force 
T=Ti-To //C 


//Rate of heat loss 

Q=T/TR 

printf(”Rate of heat loss is %fW\n”,Q) 
//(b)thermal conductivities of insulating layer 

k=til/(A*TRi) 

printf (” thermal conductivities of insulating layer 
is %f W/m C” ,k) 


Scilab code Exa 2.3 fraction of resistance 


//Example 2.3 
//Page no. 19 
// Given 

//Length & Inside rdius of gas duct 
L = 1 //nr 

r i = 0.5 / /:nr 

//Properties of inner and outer layer 
ki = 1.3 //W/m C, thermal conductivity 


bricks 

ti=0.27 
ko =0.92 
bricks 
to =0. 14 
Ti=400 
To = 65 


//nr, inner layer thickness 
//W/nr C, thermal conductivity 

//nr, outer layer thickness 
//C, inner layer temp. 

//C, outer layer temp. 


of inner 


of special 
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/ / calculation 

r_ = ri + ti //m, outer radius of fireclay brick 

layer 

ro = r_ + to I/ m, outer radius of special brick layer 

//Heat transfer resistance 

//Heat transfer resistance of fireclay brick 

Rl = (log(r_/ri))/(2*°/„pi*L*ki) 

//Heat transfer resistance of special brick 

R2 = (log(ro/r_) ) /(2*°/„pi*L*ko) 

//Total resistance 

R=R1+R2 

//Driving force 

T=Ti-To 

//Rate of heat loss 

Q=T/(R) 

printf(”Rate of heat loss is %fW’,Q) 

//interface temp. 

Tif=Ti-(Q*Rl) 

printf (” interface temp, is %f C”,Tif) 

//Fractional resistance offered by the special 
brick layer 

FR=R2/(R1+R2) 

printf (” Fractional resistance offered by the 
special brick layer is %f ”,FR) 


Scilab code Exa 2.4 Calculate Temprature 


//Example 2.4 

//Calculate (a) hot end temprature’ 

//(b) temprature fradient at both the ends 

//(c) the temprature at 0.15m away from the cold end 

//Given 
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dl=0.06 

//m, one end diameter of steel 

rod 

d2 =0. 12 

/ /m,other 

end diameter of steel 

rod 

1=0.2 

/ /m length 

of rod 


T2 =30 

// C, temp . 

at end 2 


Q =50 

//W, heat 

loss 


k = 15 

/ /W/m c , 

thermal conductivity 

of r 


//NUMERIC PART 

//T = 265.8 -(7.07/(0.06 -0.15*x) ).(a) 

//(a) 

xl=0 

// from eq . ( a) 

T1=265.8-(7.07/(0.06-0.15*xl)) 
printf(”The hot end temp, is %f C\n”,Tl) 

//(b) from eq . ( i ) 

C=50 //integration constant 

// from eq . ( i ) 

D1 = -C/(°/„pi *dl ~2*k) //D=dT/dx, tempratnre gradient 

printf(”The temprature gradient at hot end is %f C/m 
\n”,D1) 

/ / similarly 

D2 = -1179 //at x= 0.2m 

printf(”The temprature gradient at cold end is %f C/ 

m\n”,D2) 


//(c) 

x2=0.15 //m, given, 

x3 = l-x2 //m, section away from the cold 

end 

// from eq . ( a) 

T2=265.8-(7.07/(0.06-0.15*x3)) 

printf(”the temprature at 0.15m away from the cold 
end is %f C”,T2) 


Scilab code Exa 2.5 calculate refrigeration requirement 
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// Exaple2 . 5 
//Page no.24 
//Given 


// inside and 

outside 

diameter and Temp. of 

sphorical 

vessel 


do = 16 



t =0 . 1 



Ri=do/2 

//m, 

inside radius 

Ro=Ri+t 

//m. 

outside radius 

To =27 

//c, 


n 

•H 

H 

// C 


k=0.02 

/ /W/m 

. C, thermal conductivity 

layer 




//from eq. 2.23 the rate of heat transfer 

Q=(Ti-To)*(4*%pi*k*Ro*Ri)/(Ro-Ri) 
printf(”the rate of heat transfer is %f W\n’' ,[j) 
//Refrigeration capacity (RC) 

// 3516 Watt= 1 ton 
RC=-Q/3516 

printf (” Refrigeration capacity is %f tons”,RC) 


Scilab code Exa 2.6 calculate temp gradient 
//Example 2.6 

//Calculate the temprature gradient at each end of 
the rod 

//and the temprature midway in the rod at steady 
state 
//Given 

d = 0.05 //m, diameter of rod 

1=0.5 //m, length of rod 

Tl=30 //CTemp. at one end (1) 

T2=300 //C, temp at other end (2) 

T = poly (0 , ’T’) 

k = 202 + 0.0545*T //W/mC thermal conductivity of 
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//CALCULATION OF HEAT FLUX 
x 1 = 1/2 //m, at mid plane 

//temprature distribution , 

//comparing with quadratic eq. ax~2+bx+c 
//and its solution as x=(—b+sqrt(b"2 — 4*a*c))/2*a 
a = l.35*10 ~-4 
b = l 

c=-(564*xl+30.1) 

T=(-b+sqrt (b~2-4*a*c))/(2*a) 

printf(”the temprature midway in the rod at steady 
state is %f C\n”,T) 

//Temprature gradient at the ends of the rod 
x2=0 //m, at one end 

al=1.35*10~-4 
bl = 1 

cl=-(564*x2+30.1) 

Tl = (- bl + sqrt (bl~2-4*al*c1))/(2*al) 
kl=202 + 0.0545* T1 

01 = 113930 //integration constant from eq . 

(1) 

TGl = Cl/kl //C/W, temprature gradient , dT/ 

dx 

//similarly 

x3 = 0.5 

a2 = 1.35*10~-4 
b2 = l 

c2=-(564*x3+30.1) 

T2 = (-b2 + sqrt (b2~2-4*a2*c2))/(2*a2) 

k2 = 202 + 0.0545* T2 

TG2=Cl/k2 

printf (” Temprature gradient at one end of the rod is 
%f C/W\n”,TGl) 

printf (” Temprature gradient at other end of the rod 
is %f C/W’ , TG2) 
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Scilab code Exa 2.7 surface emp and maximun temp 
//Example 2.7 

//(a) what are the surface tempratures and average 
temp . of wall . 

//(b) calculate the maximum temp, in the wall and its 
location 

//(c) calculate the heat flux at the surface. 

//(d) if there is heat generation then what is the 
// average volumetric rate of heat generation? 
//Given 

x = poly (0 , ’x ’) 

/ / tem prat u re distribution in wall 

T=600+2500*x-12000*x~2 

t=0.3 //m, thickness of wall 

k = 23.5 //W/m c thermal conductivity of 

wall 

/ / Calculation 

x 1 =0 

Tl=600+2500*xl-12000*xl~2 //C, at surface 

x2 =0.3 

T2=600+2500*x2-12000*x2~2 //C, at x=0.3 

Tav = l/t*integrate( ’600 + 2500*x — 12000*x~2 ’ , ’x’ ,0,0.3) 
printf(”At the surface x=0, the temp, is %f C\n”,Tl) 
printf(”At the surface x = 0.3m, the temp, is %f C\n” , 
T2) 

printf(”Rhe average temprature of the wall is %f C” , 
Tav) 


//(b) 

D=derivat (T) //D=dT/dx 

//for maximum temprature D=0 
x3=2500/24000 
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printf(”The maximum temprature occurs at %f m\n”,x3) 

Tmax = 600 + 2500* x3-12000* x3 ~2 

printf (’’The maximum temp, is %f C\n”,Tmax) 


//(c) 

Dl=2500-24000*xl 

gradient 
Hfl=-k*Dl 

surface (x = 0) 
D2=2500-24000*x2 
gradient 
Hf2=-k*D2 

surface (x = 0.3) 
printf (’’ heat flux at 
) 

printf (” heat flux at 
Hf 2) 


//at x = 0, temprature 
//W/m"2, heat flux at left 
//at x = 0.3, temprature 
//W/m"2, heat flux at right 
left surface is %f W/m“2\n” ,Hf 1 
right surface is %fW/m~2\n”, 


//(d) 

Qt = Hf2-Hfl //W/m" 2 , total rate of heat 

loss 

Vw = 0.3 //m"3/m A 2, volume of wall 

per unit surface area 

Hav=Qt/Vw //W/m~3, average volumetric 

rate 

printf (’'The average volumetric rate if heat 
genaration is %fW/nT3 ”,Hav) 


Scilab code Exa 2.8 percentage of total heat 
//Example 2.8 

//Derive equtations for temprature distribution, 
//calculate the maximum temp, in the assembly 
//Given 

ka = 24 //W/mC thermal conductivitiy of 
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material A 


t A = 0. 1 
kB=230 
kC=200 
tB = 0.1 
tC=0.1 
TBo = 100 
TCo = 100 
Q = 2.5 * 10 ~ 5 


//m, thickness of A material 
//W/niC thermal conductivity of metl B 
//W/niC thermal conductivity of metal C 
//m, thickness of B metal 
//m, thickness of C metal 
//C, outer surface temp, of B wall 
//C, outer surface temp, of C wall 
//W/nU3, heat generated 


//NUMERIC PART 

//Temprature distribution in A, B and C 

x = poly (0 , ’x ’) 

TA=-5208*x~2+2175*x-74.5 


TB=100+96.6*x 


TC=155.2-14*x 


//position of maximum temprature x, 

D=derivat (TA) 

//At D=0 
x=2175/10416 

printf(”The maximum temp, will occur at a position 
%f m\n” , x) 
xl = x 

TA=-5208*xl~2+2175*xl-74.5 

printfC'The maximum temprature is %f C”,TA) 


Scilab code Exa 2.9 temprature distribution 

//Example 2.9 

//(a) derive eq. for temprature distribution 
//(b) find the maximum temp. 

//Given 

d i = 0 . 15 / / m, inner diameter 

do = 0.3 //m, outer diameter 

Ql = 100*10~3 //W/,m~ 3 , inner rate of heat generation 
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Q2=40*10~3 //W/m* 3, outer rate of heat generation 

Ti = 100 //C, temp . at inside surface 

To=200 //C, temp, at outside surface 

kl=30 //W/m C, thermal conductivity of 

material for inner layer 
k2 = 10 //W/m C, thermal conductivity of 

material for outer layer 


/ / Calculation 

//Tl = 364+100*log(r) -833.3*r~2 (1) 

//T2 = 718 + 216*log(r) —1000*r~2 (2) 

//(b) from eq . 1 

r = sqrt (100/2*833.3) 

printf (’’This radial position does not fall within 
layer l.\n Therefore no temprature maximum occurs 
in this layer.”) 

II similarly 

printf (” Similarly no temprature maximum occurs in 
layer 2.\n”) 

ro = di / / ixl , outer boundary 

Traax=To 

printf (’’The maximum temprature at the outer boundary 
is %f C”,Tmax) 
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Chapter 3 

Heat transfer coefficient 


Scilab code Exa 3.1 CALCULATE TIME REQUIRED 


//Example 3.1 

//calculate the time required for reduction . 

//Given 

//m, initial diameter of iceball 
//C, room temp. 

//ice ball temp. 

//W/rrU2 C, heat transfer coefficient 
//% for reduction 
//kg/nU3, density of ice 
//j/kg, latent heat of fusion 
3 //kg, mass of ice ball 

//rate of melting=—drn/dt 

//rate of heat adsorption =—4*%pi* r A 2* rho * dr / dt * 
lamda 

//at initial time t=0 

Cl=di/2 //constant of integration 

//if the volume of the ball is reduced by 40% of the 
original volume 
r=((l-x/100)*(di/2)~3)“(1/3) 

//time required for melting using eq. 1 
t=(di/2-r)/(h*(Tl-T2)/(rho*Lv)) 


di=0.06 

T1 =30 

T2 =0 

h=l1.4 

x=40 

rho =929 

Lv=3.35*10~5 

// m=4/3*%pi*r 
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printf(”The time required for melting the ice is %f 
s” ,t) 


Scilab code Exa 3.2 TIME FOR HEATING COIL 


//Example 3.2 

//calculate the time required for the 
//Given 


heating coil . 


P=l*10~3 
V = 220 

d=0.574*10~-3 
R=4.167 
Tr = 21 
h = 100 

coefficient 

rho =8920 
cp = 384 
percent =63 
/ / Calculation 
R_ = V ~ 2/P 

resistance 
1=R_/R 
A = 0 /„pi *d*l 
Tf=P/(h*A)+Tr 
dtf=Tf-Tr 


//W, electrical heating capacity 
//V, applied voltage 
//m, diameter of wire 
/ / ohm , electrical resistance 
//C, room temp. 

//W/m"2 C, heat transfer 

//kg/m" 3, density of wire 
//j /kg C, specific heat of wire 
//%, percent of the steady state 

/ /ohm, total electrical 


//m, length of wire 
//m"2, area of wire 
//final temp . 

//C. steady state temp, rise 
//temp, of wire after 63% rise 

T=Tr+(percent/100)*dtf 

//rate of heat accumulation on the wire 
//d/dt(m*cp*T) (1) 

// rate of heat loss 

//h*A*(T—Tr).( 2 ) 

//heat balance eq. (1)=(2) 

m = y„pi *d ~ 2* 1 * rho/4 //kg. mass of wire 

//integrating heat balance eq. 

t = integrate( ’ 1/ ((P / (m* cp ) ) — ((h*A) / (m* cp ) ) * (T—Tr) ) ’ , ’ 
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T’,21,322) 

printf (” The time required for the heating coil is %f 
s”,t) 


Scilab code Exa 3.3 Steady State temprature distribution 


//Example 3.3 

//(a) calculate the heat transfer coefficient 
//(b) what can be said about the same at the other 
surface of wall. 

//(c) what is average volumetric rate of heat 
generation 
//given 


t=0.2 

x=poly (0,’x’) 
T=250-2750*x~2 
k=l.163 

material 
Ta = 30 


//nr, thickness of wall 

//position in the wall 

//C, steady state temp, distribution 

//W/m C, thermal conductivity of 

//C, ambient temp 


/ / calculation 

//(a) at x = 0.2 let T=T1 at x=xl 

xl=0.2 

Tl=250-2750*xl~2 
//let D=dT/dx 

D=derivat (T) 

D = -5500*0.2 //C/m, at x = 0.2 

h=-k*D/(Tl-Ta) 

printf (” the heat transfer coefficient is %f W/m‘2 C 
, \n”,h) 


//(b) at other surface of wall, x=0=x2 (say) 

x2 =0 

a=-5500*0 

printf (”So there is no heat flow at other surface of 
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the wall \ n ” ) 


26 

27 //(c) 

28 A = 1 //m A 2 , area 

29 Vw=A*xl //m' 3 , volume of wall 

30 HL = h*(Tl-Ta) / /W, heat loss from unit area 

31 Vav = HL/xl 

32 printf (” average volumetric rate of heat generation 

is %fW/m'3”,Vav) 



Scilab code Exa 3.4 

THICKNESS OF INSULATION 

1 

clc ; 



2 

clear ; 



3 

//Example 3.4 



4 

//calculate the 

thickness of insulation 

5 

//and the rate of 

heat loss per meter length of pipe 

6 

//Given 



7 

id = 97*10 ~ -3 


//m, internal diameter of steam 


pipe 



8 

od = 114*10 ~ -3 


//m, outer diameter of steam pipe 

9 

i-i 

ii 

00 

o 


//bar, absolute pressure os 


saturated steam 


10 

Ti=234 


//C, temp, at 30 bar absolute 


pressure 



11 

Ts =55 


//C, skin temp. 

12 

o 

CO 

ii 

o 

H 


// C, ambient temp. 

13 

kc =0 . 1 


//W/m C, thermal conductivity of 


wool 



14 

kw=43 


//W/m C, thermal conductivity of 


pipe 



15 

h = 8 


//W/m A 2 C, external air film 


coefficient 



16 

L = 1 


//m, assume length 
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/ / Calculation 

ri = id/2 // m, 

r 1 = (114*10“-3)/2 //m,outer radius of steam 

pipe 

//thermal resistance of insulation 
// Ri=log ( ro / rl ) / (2* %pi*L*kc ) 

//Thermal resistance of pipe wall 
Rp=log (rl/ri)/(2*%pi*L*kw) 

/ /RT=Ri+Rp 

DF=Ti-Ts //C, driving force 

//At steady state the rate of heat flow through the 
insulation 

// and the outer air film are equ 
//by trial and error method : 

def f ( ’ [x] = f ( ro ) ’ , ’x=(Ti-Ts) /( log ( ro/rl ) /kc+log (rl/ri 
) /kw) — (h* ro * (Ts—To) ) ’) 

ro = fsolve (0.1 , f) 

th=ro-rl //m, required thickness of 

insulation 

Q=2*°/ 0 pi*ro*h*L*(Ts-To) 

printfC'The rate of heat loss is %fW,”,Q) 


Scilab code Exa 3.5 8 percent SOLUTION OF ALCOHOL 

//Example 3.5 
//calculate 

//(a) effective thickness of air and liquid films. 
//(b)the overall heat transfer coefficient based on 
i . d of pipe . 

//(c)the overall heat transfer coefficient based on 
od of insulation. 

//(d)the percentage of total resistance offered by 

air film . 
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//(e) the rate of heat loss per meter length of pipe. 
//(f) insulation skin temp. 


//given 

wl =8 
w2 = 92 
kl=0.155 
alcohol 
k2 =0.67 
water 
ka=0.0263 
kw=45 
wall 

ki=0.068 
id=53*10~-3 
od=60*10~-3 
t =0.04 
hi=800 
ho = 10 
L = 1 
T1 =75 
T2 =28 

/ / calculation 

//(a) 


//%, solubility of alcohol 

//%, solubility of water 

//W/m C, thermal conductivity of 

/ fW /m C thermal conductivity of 

//W/m C thermal conductivity of air 
//W/m Cthermal conductivity of pipe 

//W/m C , thermal cond . of glass 
//m, internal diameter of pipe 
//m, outer diameter of pipe 
//m, thickness of insulation 
//W/m"2 C, liquid film coefficient 
//W/m"2 C, air film coefficient 
//nr, length of pipe 
//C, initial temp . 

I/C, ambient air temp. 


km=(wl/100)*kl+(w2/100)*k2-0.72*(wl/100)*(w2/100) 
*(-(kl-k2)) 


deli=km/hi 

film 

// m > 

effective 

thickness 

of 

liquid 

delo=ka/ho 

//nr, 

effective 

thickness 

of 

air film 


pr int f (” e f f e c t i v e thickness of air is %f nmr”,deli 


*10~3) 


printf (”effective 
.” ,delo *10~3) 

//(b) 

Ai=2*%pi*id/2*L 

ri=id/2 

r_=od/2 

ro=r_+t 


thickness of liquid films is %f mm 


//nr "2, inside area 
//nr, inside radius of pipe 
//nr, outside radius of pipe 
//nr, outer radius of insulation 
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38 Ao= 2*%pi * ro*L //m~2, outer area 

39 //from eq. 3.11, overall heat transfer coefficient 

40 Ui=l/(l/hi+(Ai*log(r_/ri))/(2*%pi*L*kw)+(Ai*log(ro/ 

r_) ) /(2*°/ 0 pi*L*ki)+Ai/(Ao*ho) ) 

41 printf(”the overall heat transfer coefficient based 

on i.d of pipe is %f W/nT2 C”,Ui) 

42 

43 //(c) 

44 / / f r i m e q . 3.14 

45 Uo=Ui*Ai/Ao 

46 printf (” the overall heat transfer coefficient based 

on od of pipe is %f W/m‘2 C”,Uo) 

47 

48 //(d) 

49 R = l/(Ui*Ai) //C/W, total heat transfer 

resistance 

50 Rair = 1/ ( Ao *ho ) //C/W, heat transfer resistance 

of air film 

51 p = Rair/R 

52 printf (’’the percentage of total resistance offered 

by air film, is %f percent” ,p*100) 

53 

54 //(e) 

55 Q=Ui*Ai*(Tl-T2) 

56 printf (’’Rate of heat loss is %fW’,Q) 

57 

58 //(f) 

59 Ts=Uo*Ao*(Tl-T2)/(ho*Ao)+T2 

60 printf (” insulation skin temp, is %f C”,Ts) 


Scilab code Exa 3.6 Insulated flat headed 

1 // Example 3.6 

2 //calcu;ate the temp, of the liquid entering the 

bank . 
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//also calculate the insulation skin temp, at the 
flat 


//top surface and at the cylindrical surface. 
//Given 


id=l.5 
h = 2.5 
t1=0.006 
t2 =0.04 
Ta = 25 
T1 =80 
cp = 2000 
liquid 


//m, internal diameter of tank 
//m, height of tank 
//m, thickness of wall 
/ /m, thickness of insulation 
//C, ambient temp. 

//C, outlet temp, of liquid 
// j /kg C, specific heat of 


FR = 7 00/3600 


//KG/s , Liquid flow rate 


/ / Calculation 

ri=id/2+tl 

insulation 

ro=ri+t2 

insulation 
ki=0.05 

of insulation 


//m, inner radius of 
//m, outer radius of 
//W/m C, thermal conductivity 


hc=4 //W/m'2 C, heat transfer 

coefficient at cylindrical surface 
ht=5.5 //W/uT2 C, heat transfer 

coefficient at flat surface 
l = h + tl + t2 // m, height of the top of 

insulation 
//fromm eq . 3.10 

//heat transfer resistance of cylindrical wall 
Rc = log(ro/ri) / (2*%pi*l*ki) + l/ (2*°/„pi*ro*l*hc) 

//heat transfer resistance of flat insulated top 
surface 

Ri = (1/ (°/„pi *ro ~2) ) *( (ro-ri) / ki + l/ht) 
tdf=Tl-Ta //C, temp, driving force 

Q = tdf/Rc + tdf/Ri //W, total rate of heat loss 

Tt=Q/(FR*cp)+Tl //C, inlet temp, of liquid 

printf (” Inlet liquid temp, should be %f C \n’ : ,Tt) 

Ql = tdf/Ri //W, rate of heat loss from flat surface 
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32 T1 = Q 1 / (°/„pi *ro “2*ht) +Ta 

33 printf (” the insulation skin temp, at the flat top 

surface is %f C \n”,Tl) 

34 / / similarly 

35 T2 =38 

36 printf (” similarly the insulation skin temp at 

cylindrical surface is %f C”,T2) 


Scilab code Exa 3.7 rate of heat transfer 

1 //Example 3.7 

2 //what is the heat imput to the boiling. 


3 // Given 

4 id=2.5*10"-2 

glass tube 

5 t=0.3*10“-2 


6 1=2.5 

wire 

7 L = 0.12 

covered with heating 

8 Re=16.7 

resistance 

9 ti=2.5*10~-2 

insulation 

10 kg = 1.4 

conductivity of glass 

11 ki=0.041 

conductivity of 

12 T1 =91 

liquid 

13 T2 =27 

14 ho = 5.8 

film coefficient 

15 V = 90 

16 


//m, internal diameter of 

//nr, thickness of wall 
//nr, length of niclrronre 

//nr, length of steel 
coil 

/ / ohm , electrical 
//nr, thickness of layer of 
//W/nr C, thermal 
//W/nr C, thermal 


insulation 

//C, boiling temp, of 


//C, ambient tenrp. 
//W/nr "2 C outside air 

//V, voltage 
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17 

/ / Calculation 





18 

Rc=Re*l 


/ / ohm , resistance 

of 



heating coil 





19 

Q = V ~ 2/Rc 


//W, rate of heat 




generation 





20 

ri=id/2 


// m, inner radius 

of 

glass 


tube 





21 

r_=ri+t 


//nr, outer radius 

of 

glass 


tube 





22 

ro=r_+ti 


//nr, outer radius 

of 



insulation 





23 

//heat transfer 

resistance of glass wall 



24 

Rg=log(r_/ri)/(2* %pi*L*kg) 



25 

//combined resistance 

of insulation and outer 

air 


film 





26 

Rt = log(ro/r_) / (2*%pi*L*ki) + l/ (2*°/ 0 pi*ro*L*ho) 



27 

//Rate of heat 

input 

to the boiling liquid in 

s t e e1 = 


Ql=(Ts—Tl)/Rg 

28 //Rate of heat loss through insulation ,Q2=(Ts—To) / ( 

Rt) 

29 // Q1+Q2=Q 

30 Ts=(Q+ Tl/Rg +T2/Rt)/(1/Rg +1/Rt) 

31 Ql=(Ts-Tl)/Rg 

32 Q2 = Q-Q1 

33 printf(”the heat impnt to the boiling, is %fW’,Ql) 


Scilab code Exa 3.8 A 10 gauge electrical copper 

1 //Example 3.8 

2 //determine (a) maximum allowable current 

3 //(b) the corresponding remp. at the centre of wire 

and 

4 //at the outer surface of insulation 

5 // Given 

6 ri = 1.3*10 ~-3 


//m, radius of 10 gauge wire 
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t=l.3*10"-3 
insulation 
Ti =90 
To =30 
h = 15 

coefficient 

km=380 
copper 

kc =0.14 

rubber(insulation ) 

Rc = 0.422/100 

resistance of copper 

//NUMERIC CALCULATIONS 
T cmax = 90 

insulation 


//nr, thickness of rubber 

//C, temp. Of insulation 
//C, ambient temp. 

//W/nU2 C, air film 

//W/m C, thermal cond. of 

//W/m C, thermal cond. of 

/ /ohm/m, e 1 e t r i c a 1 
wire 

//X, the maximum temp, in 
//m, outside radius of 10 


ro=ri+t 

gauge wire 

Sv=((Tcmax-To)*(2*kc/ri~2))/ (log(ro/ri)+kc/ (h*ro) ) 
//from eq . ( xi i ) , Sv=I "2*rho/(%pi*ri "2) 

I = (°/ 0 pi*ri ~2*Sv/Rc) ~0.5 //A, Current strength 

printf (’’maximum allowable current is %f A\n”,I) 


//(b) at r=0 

Tm=To+(ri~2*Sv/2)*(l/km+(log(ro/ri))/kc+l/(h*ro)) 
printf (” rernp . at the centre of wire is %f C\n”,Tm) 

//at r=ro 

Tc=30+(ri~2*Sv/(2* kc))*(kc/(h*ro)) 
printf (’’The temprature at the outer surface of 
insulation is %f C”,Tc) 


Scilab code Exa 3.9 Heat generating slab A 
//Example 3.9 

//(a) calculate the temp, at the surface of slab A. 
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//what is the maximum Temp, in A. 

//(b) deter mine the temp, gradient at both the 
//surfaces of each of the slabs A,B and C. 
//(c)calculate the value of hi & h2. 


//Given 

tA = 0.25 
tB = 0.1 
tC=0.15 
kA = 15 

slab A 
kB = 10 

slab B 
kC = 30 

slab C 

x = poly (0 , ’x ’) 
B 


//m, thickness of slab A 

//m, thickness of slab B 

//m, thickness of slab C 

//W/m C, thermal comductivity of 

//W/m C, thermal comductivity of 

//W/m C, thermal comductivity of 

//m, distance from left surface of 


//Temprature distribution in slab A 

TA=90+4500*x-11000*x~2 


Tl=40 //C, fluid temp. 

T2=35 //C, medium temp. 


/ / calculation 

//(a) 

xl = tB 

TAl=90+4500*xl-11000*xl~2 

//similarly at the right surface 

x2=tA+tB 

TA2=90+4500*x2-11000*x2~2 
//let dTA/dx=D 
D=derivat (TA) 

D = 0 //for maximum temp. 

x3=4500/22000 

TAmax=90+4500*x3-11000*x3“2 

printf(”At x = 0.1 the temp, at the surface of slab A 
is %f C\n”,TA1) 

printf(”At x = 0.35 the temp, at the surface of slab A 
is %f C\n”, TA2) 


30 


35 printf (” the maximum Temp, in A occurs at %f m\n” , 

x3) 

36 printf (” the maximum Temp, in A is %f TArnax \n” , 

TAmax) 


37 


38 //(b) 

39 //At the interface 2 

40 Dl=4500-2*11000*xl //C/W, Dl=dTA/dx , at x = 0.1 

41 //At the interface 3 

42 D2=4500-2*11000*x2 //D12=dTA/dx , at x = 0.35 

43 //Temprature gradient in slab B and C 

44 //by using the continuity of heat flux at interface 

( 2 ) 

45 D3 = -kA *D1 / (-kB) //D3=dTB/dx, at x = 0.1 


46 //at interface (1) 

47 D4 = D3 //D4=dTB/dx at x=0 

48 //similarly 

49 D5 = -1600 //C/W, dTB/dx, x = 0.35 

50 D6 = D5 / / at interface 4 


51 

printf (” temp . 

is %f C/W\ 

gradient 
n”,D1) 

at 

interface 

2 

of 

the 

slabs 

A 

52 

printf (” temp . 
is %f C/W\ 

gradient 
n” , D2) 

at 

interface 

3 

of 

the 

slabs 

A 

53 

printf (” temp . 
is %f C/W\ 

gradient 
n”,D3) 

at 

interface 

2 

of 

the 

slabs 

B 

54 

printf (” temp . 
is %f C/W\ 

gradient 
n”,D4) 

at 

interface 

1 

of 

the 

slabs 

B 

55 

printf (” temp . 

Cis %f C/W\ 

gradient 
n” , D5) 

at 

interface 

3 

of 

the 

slabs 


56 

printf (” temp . 

gradient 

at 

interface 

4 

of 

the 

slabs 

C 


is %f C/W\n” , D6) 

57 

58 //(c) 

59 //from D3 = 3450 and TB=beetal *x+beeta2 

60 beetal=3450 

61 beeta2=85 

62 x = 0 

63 TB=beetal*x+beeta2 
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/ / s i m i 1 a r y 

TC = 877.5-1600*x 
hl=-kB*D4/(Tl-TB) 

//similarly 
h2 = 1129 

printf (” The heat transfer coefficient at one 

surface of solid fluid interface is %fW/nT2 C\n” , 
hi) 

printf (’’The heat transfer coefficient at other 

surface of solid fluid interface is %fW/m" 2 C” , h2 
) 


Scilab code Exa 3.10 percentage increase in rate 


//Example 3.10 

//calcuate the percentage increase in the rate of 
heat transfer 

//for the finned tube over the plain tube. 

// Given 

//m, actual internal dia of pipe 
//m, wall thickness 
II no . of longitudinal fins 
//m, thickness of fin 
//m, breadth of fin 
//W/m C, thermal conductivity of 


id=78*10~-3 
tw=5.5*10"-3 
nl =8 

tf=1.5*10"-3 
w=30*10"-3 
kf =45 

fin 
Tw=150 
To =28 
h = 75 

coefficient 


//C, wall temp. 

//C, ambient temp. 

//W/m"2C, surface heat transfer 


/ / Calculation 
// from eq . 3.27 

e = sqrt (2*h/(kf*tf)) 

n= ( 1/ ( e*w) ) *tanh (e*w) //efficiency of fin 


32 



19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

1 

2 

3 

4 

5 

6 

7 


//nr, length of fin 
//nT2, area of single fin 
//m A 2 total area of fin 

//W, maximum rate of heat 

//W, actual rate of heat 

/ / m A 2 , area of contact of fin 


L = 1 

Af =2*L*w 
Atf=nl*Af 

Qmax=h*Atf*(Tw-To) 
transfer 
Qa = n* Qmax 
transfer 
Afw=L*tf 

with pipe wall 

Atfw=Afw*nl //uT2 , area of contact of all 

fin with pipe wall 

ro = id/2 + tw // m, outer pipe radius 

A = 2* °/ 0 pi *L*ro //m"2 area per meter 

Afree=A-Atfw //uC2, free outside area of 

finned pipe 

//Rate of heat transfer from free area of pipe wall 
Ql=h*Afree*(Tw-To) / /W, 

//total rate of hewat gtransfer from finned pipe 

Qtotal=Qa+Ql //W 

//Rate of heat transfer fromm unfinned pipe 

Q2=h*A*(Tw-To) 
per=(Qtotal-Q2)/Q2 

printf(”the percentage increase in the rate of heat 
transfer is %f percent ”,per*100) 


Scilab code Exa 3.11 Pre stresed multilayered shell 

//Example 3.11 
/ / Calculate 

//(a) Is there any thermal contact resistance at the 
interface between the layer? 

//(b) if so calculate the contact resistance and 
//express it in contact heat transfer coefficient 
//(c) Calculate the temp. jump. 
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8 

//Given 




9 

id=90*10~-2 

// m , 

internal diameter 

of steel 

10 

od=110*10~-2 

// m > 

outer diameter of 

steel 

11 

Ti=180 

//c, 

inside temp, of steel 

12 

To = 170 

//C, 

outside temp. of 

steel 

13 

h- 

co 

ii 

//W/m 

C, thermal conductivity of 


alloy 




14 

15 

16 

Q = 5.18*10 ~ 3 

// w, 

Rate of heat loss 


/ / calculation 




17 

ri=id/2 

// m , 

inside radius of 

shell 

18 

ro=od/2 

// m > 

outside radius of 

shell 

19 

r_ =0.5 

// m , 

boundary between 

the layers 

20 

L = 1 

// m > 

length of shell 



21 //Rate of heat transfer in the absence of contact 

resistance 

22 Ql=2*%pi*L*k*(Ti-To)/ (log(ro/ri) ) 

23 printf(”Rate of heat transfer in the absence of 

contact resistance is %f KW\n” , Q1 /1000) 

24 printf(”The actual rate of heat loss is 5.18kW is 

much less than this value.\n So there is a 
thermal contact resistance at the interface 
between the layers \n” ) 

25 


26 //(b) 

27 Ri = ( log (r_/ri ) / (2* °/„pi *L*k) ) //C/W, Resistance of 

inner layer 

28 Ro = ( log (ro/r_ ) / ( 2* °/ 0 pi *L*k) ) //C/W, Resistance of 

outer layer 

29 Rc = ((Ti-To ) / (Q) ) - (Ri + Ro) //C/W, contact 

resistance 

30 printf(”The contact resistance is %f C/W\n”,Rc) 

31 Ac =2* %pi *L* r_ //uT2, area of contact 

surface of shell 

32 hc = l/(Ac*Rc) / /W/m' 2 c, contact heat 

transfer coefficient 

33 printf (” contact heat transfer coefficient is %f W/m 

A 2 C \n”,he) 
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//(c) 

dt=Q/(hc*Ac) 

printf (’’The temprature jump is %f C” , dt ) 


Scilab code Exa 3.12 critical insulation thickness 


// Example 3.12 
// calculate the 

d = 5.2 * 10 ~-3 
ri=d/2 
kc =0.43 
PVC 
Tw = 60 
h=l1.35 
To =21 

/ / calculation 

Ro=kc/h 

insulation 


critical thickness . 

//m, diameter of copper wire 
//inner radius of insulation 
//W/m C, thermal conductivity of 

//C, temp . Of wire 

//W/m"2 C, film coefficient 

//C, ambient temp. 

//m, critical outer radius of 


t=Ro-ri 

printf (’’the critical thickness is %f mm” , t * 10~3) 


Scilab code Exa 3.13 critical insulation thickness 


//Example 3.13 
// calculate the 

d=15*10~-2 
t = 10 * 10"-2 
ki=0.035 

insulation 
h = 10 

coefficient 
/ / calculation 


critical insulation thickness, 
//m, length of steam main 
/ / m, thickness of insulation 
//W/m C, thermal conductivity of 

//W/m"2 C, heat transfer 


35 




8 

9 

10 

11 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 


/ / from eq . 3.2 9 

ro=ki/h 

printf (”ro= %f cm \n” ,ro* 10~3) 

printf (” Radius of bare pipe is larger than outer 
radius of insulation \n So critical insulation 
thickness does not exist ”) 


Scilab code Exa 3.14 optimum thickness 


//Example 3.14 
//calculate the 
//Given 

Ti=172 
To =20 
Cs =700 
Lv=487 
ho=10.32 

coefficient 
kc =0.031 

insulation 


optimum thickness . 

//C, saturation temp. 

//C, ambient temp. 

//per ton, cost of steam 
//kcal/kg, latent heat of steam 
//kcal/h m"2 C, outer heat transfer 

//W/m C, thermal conductivity of 


n=5 //yr, service life of insulation 

i=0.18 //Re/(yr)(Re), interest rate 

/ / Calculation 

di=0.168 //m, inner diameter of insulation 

//Cost of insulation 

Ci=17360-(1.91*10~4)* di //Rs/m"3 

Ch=Cs/(1000*Lv) //Rs/cal, cost of 

heat energy in steam 

sm = l/ (l + i)+l/(1 + i)~2+l/(l + i)~ 3 + 1/(1 + i)“4+1/(1 + i) ~n 

/ / from eq . 3.3 3 

ri = di/2 / /m inner radius of insulation 

L = 1 //m, length of pipe 

/ / Pt=Ch*sm*2*%pi* ri*L*( l/(((ri/kc)*(’log(ro/ri) ’)) + 
ri / ( ho* ro ) ) ) *7.2*10"3*( Ti—To)+%pi * ( ro "2 — r i "2)*L* 
Ci 
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22 //On differentiating , dpt / dro = — 957.7* ((1/ro ) 

- (0.003/ro "2) ) /(log (ro) +(0.003/ro) +2.477) ~2 

23 deff ( ’ [x]=f (ro) \ ’ x= -957.7* ((1 / ro ) - (0.003/ro * 2) ) / ( 

log (ro) +(0.003/ro) +2.477) ~2+98960* ro ’) 

24 ro = fsolve(0.1,f) 

25 t = ro-ri 

26 printf(”The optimum insulation thickness is %f mm” , t 

*1000) 
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Chapter 4 

Forced Convection 


Scilab code Exa 4.2 Air flow over a flat plate 


//Example 4.2 
//Determine 

//(a) local heat transfer coefficient 
//(b)the average heat transfer coefficient 
//the rate of heat loss from the surface. 


//Given 

1 = 2 

T1=150 
p = l 
T2 =30 
V = 12 


//m, length of flat surface 
//C, surface temp. 

//atm , pressure 
//C, bulk air temp. 

//m/s, air velocity 


/ / Calculation 

Tf=(Tl+T2)/2 
mu=2.131*10"-5 
k=0.031 

conductivity 
rho =0.962 
cp = 1.01 
of air 


//C, mean air film temp. 
//m~2/s, viscosity 
//W/m C, thermal 

//kg/m~3, density of air 
//kj/kg C, specific heat 
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Pr = cp* 10"3*mu/k //Prandtl no. 

Remax = l*V*rho/mu //maximum Reynold no. 

Re=5*10"5 //Reynold no. during 

transition to turbulent flow 
L_ = ( Re *mu) / (V* rho ) //m, distance from the 

leading edge 

//for laminar flow heat transfer coefficient h, 
//hl6.707*x A —(1/2) 

//(a) 

//h2 = 31.4*x"(-l/5) 

//b 

hav=22.2 

//c 

Q=hav*l*p*(Tl-T2) 

printf(”The rate of heat loss is %fW” ,Q) 


Scilab code Exa 4.3 temprature of wire 


//Example 4.3 

//what will be the temp. 

//Given 

d=7.24*10"-4 
1 = 1 
1=8.3 
R=2.625 

resistance 

V = 10 
Tb = 27 

//the properties at bulk 
mu=l.983*10"-5 
k=0.02624 

conductivity 

rho = 1.1774 

air 


of the wire at steady state 

/ /m, diameter of wire 
//m, length of wire 
//A, current in a wire 
/ / ohm /m, electrical 

//m/s, air velocity 
//C, bulk air temp. 

temp . 

//nr" 2/s, viscosity 
//W/m C, thermal 

//kg/m"3, density of 
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//kj/kg C, specific 


cp = l .0057 

heat of air 


/ / calculation 

Pr = cp * 10~3*mu/k //Prandtl no. 

Re=d*V*rho/mu // Reynold no. 

//from eq. 4.19, nusslet no. 

Nu = 0.3+(0.62*Re~(l/2)*Pr~(l/3)/(l+(0.4/Pr) "(2/3)) 
~(1/4))*(l+(Re/(2.82*10“5))~(5/8))~(4/5) 
hav=Nu*k/d //W/m"2 C, average 

heat transfer coefficient 


Q=I~2*R 

electrical heat generation 

A = °/ 0 pi *d*l 
dt=Q/(hav*A) 

T=dt+Tb 


//W, rate of 


//C,temp. difference 
//C, steady state temp 


printf(”The steady state temprature is %f C\n”,T) 
//REVISED CALCULATION 

Tm=(T + Tb)/2 //C, mean air film 

temp . 

//the properties at Tin temp. 

mul=2.30*10~-5 
kl=0.0338 

conductivity 
rho1=0.878 

air 

cpl=l.014 

heat of air 
Rel=d*V*rhol/mul 
Prl=(1.014*10~3*2.30*10~-5)/kl 
no . 

//from eq. 4.19, nusslet no. 

Nul=0.3+(0.62*Rel~(1/2)*Prl~(1/3)/(l+(0.4/Prl)“ (2/3) 
)~(1/4))*(l+(Rel/(2.82*10~5))“(5/8))~(4/5) 
havl=Nul*kl/d //W/m"2 C, average 

heat transfer coefficient 

dt 1 = Q / (havl * A ) //C,temp. difference 


//nV2/s, viscosity 
/ fW /m C, thermal 

//kg/m" 3, density of 

//kj/kg C, specific 

// Reynold no . 

// Prandtl 


40 


//C, steady state 


40 T1 = dt 1 +Tb 

temp . 

41 printf ( r The recalculated value is almost equal to 

previous one . ” ) 


Scilab code Exa 4.4 Calculate Required time 

1 //Example 4.4 

2 / / Calculate 

3 //(a) what is initial rate of melting of ice. 

4 //(b)how much time would be needed to melt away 50 % 


of ice 


5 // Given 

6 di=0.04 


/ /m, diameter of ice 


ball 

7 V = 2 

8 T1=25 

9 T2 =0 

10 //the properties of air 

11 mu=l.69*10~-5 

12 k = 0.026 


//m/s, air velocity 
//C, steam temp. 


//kg/ms, viscosity 
//W/m C, thermal 


conductivity 


14 cp = 1.005 

heat 


13 rho = 1.248 


//kg/m"3, density 
//kj/kg C, specific 


15 //proper! ice of ice 


16 lamda = 334 

fusion 


//kj/kg , heat of 


17 rhoice=920 

ice 


//kg/nT3 density of 


18 

19 II calculation 

20 Pr = cp*10~3*mu/k 

21 Re=di*V*rho/mu 


/ / Prandtl no . 

// Reynold no . 


22 //from eq. 4.19, nuss let no. 
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Nu=2+(0.4*Re"0.5+0.06*Re~(2/3))*Pr"0.4 
hav=Nu*k/di //W/nT2 C, average 

heat transfer coefficient 
Ai = °/ 0 pi*di~2 //initial area of 

sphere 

Qi = Ai*hav*(T1-T2 ) //W=J / s , initial rate 

of heat transfer 

Ri=Qi/larada //initial rate of 

melting of ice 

printf (” i n i t i a 1 rate of melting of ice is %f g/s\n” 

, Ri) 

//(b) 

//mass of ice ball 4/3*%pi*r"3 
//Rate of melting= Rtm= — d/dt(m) 

//Rate of heat input required =— lamda*Rate of 
melting 

//heat balance equation 
// — larnda * (Rm)=h*4* %pi* r " 2* dt 
//integrating and solving 

rf=((di/2)"3/2)"(1/3) 

//solving eq. 3 
t1=1.355*10"-4/(8.136*10"-8) 

printf (’’The required time is is %f s \n” , round ( 1 1 ) ) 


Scilab code Exa 4.5 average time 


//Example 4.5 
//calculate the average 
//Given 

Vo=0.5 
T1=800 
T2 = 550 
Tam=500 
P = 1.2 


time of contact . 

//m/s air velocity 
//C, initial temp . 
//C, final temp . 
//C, air mean temp 
//atm , pressure 
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//the properties of solid 

dp = 0.65*10"-3 

diameter 

cps =0.196 

heat 

rhos =2550 

//Properties of air 

mu=3.6*10"-5 
k=0.05 

conductivity 

rho =0.545 

air 

cp = 0.263 

heat of air 


particles . 

//m, average particle 

//kcal/kg C, specific 

//kg/m" 3, density 

//kg/ms, viscosity 
//kcal/hm C, thermal 

//kg/m"3, density of 

//kcal/kg C, specific 


/ / calculation 

Pr=cp*mu*3600/k //Prandtl no. 

Redp=dp*Vo*rho/rau // Reynold no. 

//from eq. 4.29(b) heat transfer coefficient 

h=(k/dp)*(2 + 0.6*(Redp)"(1/2)*(Pr)" (1/3)) 

Tg=500 //C, gas temp. 

//from heat balance equation 
// — (dTs/dt)=6h/(dp*rhos*cps)*(Ts—Tg) 
t=(dp*rhos*cps/(6*h))*integrate( ’(1/(Ts-Tg))’,’Ts’ 
,550,800) 

printf(”the required contact time is %f s”,t*3600) 


Scilab code Exa 4.6 Overall heat transfer coefficient 
//Example 4.6 

//Calculate the required rate of flow of water, 
//calculate the overall heat transfer coefficient 
//Given 

rao_ = 1000 //kg/h , cooling rate of oil 

cpo=2.05 //kj/kg C, specific heat of oil 
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T1 =70 
T2 =40 
cpw=4.17 
T3 =42 
T4 = 28 
A = 3 

/ / Calculation , 


//C, initial temp, of oil 
//C, temp, of oil after cooling 
//kj/kg C, specific heat of water 
//C, initial temp, of water 
/ / C, temp, of oil after cooling 
//m"2, heat exchange area 
rate of flow of water 


mw_=mo_*cpo*(Tl-T2)/(cpw*(T3-T4)) 

printf(”the required rate of flow of water is %f kg/ 
h \n”,mw_) 

Q=mo_*cpo*( T1 -T2)/3600 //kw, heat duty 

dtl=Tl-T3 //C, hot end temp, difference 

dt2 = T2-T4 //C, cold end temp, difference 

LMTD=(dtl-dt2)/(log(dtl/dt2)) //log mean temp. 

difference 

dtm=LMTD 


U=Q*10~3/(A*dtm) 

printf(”the overall heat transfer coefficient is %f 
W/rn' 2 C”,U) 


Scilab code Exa 4.7 inlet and outlet temprature 


//Example 4.7 
// calculatemthe 
//Given 

Q =38700 
W=2000 
cp = 0.239 
nitrogen 
A = 10 
U = 70 


inlet and outlet temp, of gas. 

//kcal/h, heat duty 
//kg/h gas flow rate 
//kcal/kg C, specific heat of 

//m"2 ,heat exchanger area 

//kcal/hm"2 C, overall heat 


transfer coefficient 
n = 0.63 //fin efficiency 


/ / Calculation 
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dt = Q/(W*cp) //C, temp, difference 

//To—Ti=dt.( i ) 


dtm=Q/(U*A*n) 

/ / (To—Ti) / (log((160 — Ti)/(160—To))) = 87.8.(2) 

//solving 1 and 2 

deff ( ’ [x]= f (To) ’ , ’ x=(To—(To—dt) ) / ( log ((160 - (To-dt) ) 

/(160 —To) ) ) -87.8 ’) 

To = fsolve (100 ,f) 

Ti=To-dt 

printf(”The inlet temprature is Ti=%f C\n”, round (Ti 

)) 

printf(”The outlet temprature is To=%f C\n ,: , round (To 
)) 


Scilab code Exa 4.8 drop in temprature 


//Example 4.8 
//Calculate the 
//Given 

V = 1.8 
TI=110 
1 = 15 
t=0.02 
kc =0.12 

of insulating 
ho = 10 

coefficient 

T2 =20 

//the properties 

mu=2.55*10“-4 
k=0.685 
rho =950 
cp =4.23 

air 

di=0.035 


drop in temp, of the water. 

//m/s, velocity of hot water 
//C, initial temp . 

//m, length of pipe 

//m, thickness of insulation 

//W/mC, thermal conductivity 


layer 


/ /Win" 2 C, outside film 


//C, ambient temp, 
of water at 110 C 

//m"2/s, viscosity 
//W/m C, thermal conductivity 
//kg/m"3, density of air 
//kj/kg C, specific heat of 

//m, actual internal dia . of 
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pipe 

ri=di/2 
t1=0.0036 

1 — 1/ 4 schedule 

ro=ri+tl 

r_=ro+t 


/ / nr, internal radius 
//nr, actual thickness of 
40 pipe 

//nr, outer radius of pipe 
//nr, outer radius of 


insulation 

kw=43 //W/mC, thermal conductivity 

of steel 
/ / calculation 

Pr = cp* 10~3*mu/k //Prandtl no. 

Re=di*V*rho/mu // Reynold no. 

//from eq. 4.9, Nusslet no. 

Nu = 0.023*(Re)~0.88*Pr~0.3 


hi = Nu*k/di //W/nU2 C, average heat 

transfer coefficient 

//the overall coefficient inside area basis Ui 

Ui=l/(l/hi+(ri*log(ro/ri))/kw+(ri*log(r_/ro))/kc+ri 
/(r_ *ho)) 

Ai = °/„pi *di *1 //m~2, inside area basis 

W = 0 / 0 pi * r i ~ 2* V* rho //kg/s, water flow rate 

//from the relation b/w LMTD and rate of heat loss 
// def f ( ’ [ x] = f (To) ’ , ’ xAV* cp*lCU3* (Tl—To)—Ui* Ai *((Tl— 
To) / log ((Tl—T2) / (To—T2) ) ) ’) 

//To=fsolve(1 , f ) 


deff ( ’ [x]=f (To) ’ , ’x=(W*cp *10'3) /( Ui* Ai ) * (Tl-To) - ((Tl 
—To) / log ((Tl—T2) / (To—T2) ) ) ’) 

To = fsolve (100 ,f) 

printf(”The outlet eater temp, is %f C”,To) 


Scilab code Exa 4.9 find the tenrprature 
//Example 4.9 

//at what tenrp. does the water leave the pipe. 
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//Given 

T1 =28 
T2 = 250 
V = 10 
1=20 

mw=1*3600 
di = 4 . 1 * 10"-2 
Tm=(Tl+T2)/2 
ro = 0.0484 
//properties of 
mu=8.6*10"-4 
kw=0.528 

conductivity 
kw_ =0.528*1. 162 
rho =996 
cp = 1*4.18 
cp_ = 1 

//properties of 

mul=2.33*10"-5 
ka=0.0292 

conductivity 
rho1=0.891 
cpl=0.243 
Pr = 0.69 


//C, inlet temp. 

//C, bulk temp . 

//m/s , gas velocity 
//m, length of pipe 
//kg/h , water flow rate 
//m, inlet diameter 
//C, mean temp. 

/ / m, outside radius 
water 

//kg/ms, viscosity 
// kcal/h m C, thermal 

/[ W/ m C, thermal conductivity 
//kg/m" 3, density of air 
//kj/kg C, specific heat of air 
//kcal/kg C 
flue gas 

/ / kg /ms , viscosity 
//kcal/h m C, thermal 

//kg/m"3, density of air 
// kcal/kg C, specific heat of air 


II calculation 

A = °/„pi/4* di "2 //m"2, cross section of pipe 

Vw = l/(rho*A) //m/s, velocity of warer 

Re=di*Vw*rho/mu // Reynold no. 

Prl=cp*10"3*mu/kw_ //Prandtl no. for water 
Nu = 0.023*Re "0.8*Prl "0.4 //Nusslet no. 

//water side heat transfer coefficient hi 
hi=206* kw/di 

//gas side heat transfer coefficient ho 
a=41 //mm, i.d. schedule 

Tw = 3.7 //mm, wall thickness 

do = a + 2*Tw //mm, outer diameter of pipe 

Re 1 = do* 10"-3*V*rho1/mul // Reynold no 
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//from eq. 4.19, nusslet no. 

Nul=0.3+(0.62*Rel~(l/2)*Pr~(l/3)/(l+(0.4/Pr)~(2/3)) 

~(1/4))*(l+(Rel/(2.82*10~5))~(5/8))~(4/5) 
ho = ( Nul *ka/do ) * 10 ~ 3 //kcal/li rn' 2 C 

Uo=l/(ro/(di/2*hi)+l/ho) //kcal/h m'2 C, overall 
heat transfer coefficient 

//Heat balance 

Al = y„pi*ro*i //m62, outside area of pipe 

//from the formula of LMTD 

deff ( ’ [x]=f (T2_) ’ , ’ x^ikv* cp _ * (T2_—T1)—Uo*Al * (( T2_—T1) 
/log ((T2—T1) / (T2—T2_) ) ) ’) 

T2_ = fsolve (1 , f ) 

printf (” The exit water temp is %f K” , round ( T2_ ) ) 


Scilab code Exa 4.10 length of heat exchanger 


//Example 4.10 
//calculate the 
// Given 

dti=0.0212 
dto =0.0254 
dpi=0.035 
rao_ =500 
To2 = 110 
To 1=70 
Tw2 =40 
Twl=29 

//properties of 

cpo =0.478 
ko =0.12 
rho =850 

//properties of 

kw=0.542 

kw_ = (kw * 1.162) 


length of heat exchanger. 

//nr inner tube 
//cm , outer tube 
//cm , outer pipe 
//kh/li, cooling rate of oil 
//C, initial temo. of oil 
//C, temp, after cooling of oil 
//C, inlet temp, of water 
//C, outlet temp, of water 
o i 1 

//kcal/kg C 

//kcal/h nr C, thermal conductivity 
//kg/nr~3, density of oil 
water 

//kcal/h nr C, thermal conductivity 
//kj/kg C 
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rauw=7.1*10“-4 
cpw = l 

cpw_=cpw*4.17 
rhow=1000 

/ / calculation 

HL=mo_*cpo*(To2-Tol) 

exchanger 

raw_ = HL/(cpw*(Tw2-Twl) ) 
mw_l=mw_/(3600*10“3) 

A1 = (°/„p i /4) *(dti) “2 
Vw=mw_1/A1 
Rew=dti*Vw*rhow/muw 


//kg/ms, viscosity of water 
//kcal/kg C 
//kcal/kg C 
//kg/m“3, density 


//kcal/li, heat load of 

//kg/h water flow rate 
//nT3/s water flow rate 
//m“2, flow area of tube 
//m/s water velocity 
//Reynold no. 

//Prandtl no. 


Prw = cpw_*10“3* muw/kw_ 

Nuw = 0.023* Rew “ 0.8* Prw “ 0.4 //nusslet no. 
//water side heat transfer coefficient hi 
hi=Nuw*kw/dti 


//oil side heat transfer 

A2 = 7„pi/4* (dpi “2-dto “2) 

annulus 

Vo=mo_/(3600*rho*A2) 
de=(dpi~2-dto“2)/dto 

annulus 

Tmo=(To2+Tol)/2 

muoi1 = exp ((5550/(Tmo+273)) -19) 

of oil 

Reo=de*Vo*rho/muoil 

Pro=cpo*muoil*3600/ko //prandtl no. 


coefficient 
//m“2, flow area of 

//m/s velocity of oil 
//m, equivalent dia of 

//C,mean oil temp. 

//kg/ms, viscosity 


for oil 


//assume (1st approximation) 

Nuo =3.66 

ho=Nuo*ko/de //kcal/h m'2 c 

L = 1 //assume length of tube 

Ai = °/„pi *dti *L 
Ao = °/„pi *dto *L 

//overall heat transfer coefficient 1st 
approximation 

Uo=l/(l/ho+Ao/(Ai*hi)) 
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LMTD = ((To2-Tw2)-(To 1 -Twl))/( log ((To2-Tw2)/(To 1 -Tw1)) 

) 

Ao 1 = HL/ (Uo*LMTD) //irT 2 , heat transfer area 

Lt = Ao 1 / (°/ 0 pi * dt o ) //m, tube length 

// from eq . 4.8 

Nuo 1 = 1.86* ( Reo * Pro / (Lt/de ) ) " (1 /3) //Nusslet no. 

hol=Nuol*ko/de 

Tmw = (Twl+ Tw2)/2 //C, mean water temp, 

//balancing heat transfer rate of oil and water 

//average wall temp. Twall 

Twall=((hi*dti*(-Tmw))-(hol*dto*Tmo))/(-65.71216) 

//viscosity of oil at this temp. 

muwal1 = exp ((5550/(Twall+273))-19) //kg/ ms , 
viscosity of oil 
//Nusslet no. 

Nuo2=1.86*(Reo*Pro/(Lt/de))~(l/3)*(muoil/muwal1) 
"0.14 

ho2=Nuo2*ko/de 

Uo2=l/((l/ho2)+(Ao/(Ai*hi))) 

Ao2=HL/(Uo2*LMTD) 

Lt_=Ao2/(%pi*dto) 

printf(”The tube length is %f m”,Lt_) 


Scilab code Exa 4.11 rate of heat transfer 
//Example 4.11 

//calculate the rate of heat transfer to water. 
//Given 

Ti=260 //C, initial temp. 

Ts=70 //C, skin temp. 

St=0.15 //m, space between tubes in 

equilateral triangular arrangement 
Sd = St //space between tubes 

mu = 4.43* 10" -5 //m"2/s , momentum 
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diffusity 
k=0.0375 

conductivity 
rho =0.73 

air 

cp = 0.248 

heat of air 
V = 16 
d=0.06 

of tube 
Nt = 15 

transverse row 

N1 = 14 

longitudinal row 

N=Nl*Nt 
L = 1 

/ / Calculation 

SI = ( sqrt (3) /2) *St 
Pr=cp*mu*3600*rho/k 

air 

Pr = 0.62 
Prw=0.70 

wall temp. 70 C 
// from eq . 4.25 
Vmax=(St/(St-d))*V 
// from eq . 4.26 
Vmaxl=(St/(2*(St-d)))*V 
Redmax=d*Vraax/mu 
p=St/Sl 

p <2 

// from table 4.3 


/ fW /m C, thermal 

//kg/m" 3, density of 

//kj/kg C, specific 

//m/s , ve1o ci i t y 
//m, outside diameter 

//no. of tubes in 

//no. of tubes in 

//total no. of tubes 
//m, length 

//Prandtl no. of bulk 

//Prandtl no. of air at 


//pitch ratio 


ra = 0.6 

C = 0.35*(St/SI)~0.2 

h=(k/d)*C*(36163)"m*(Pr)“(0.36)*(Pr/Prw) "(0.25) 

// from eq . 4.28 

dt = 190*exp(-°/ 0 pi*d*N*h/(rho*V*3600*Nt*St*cp)) 

LMTD=((Ti-Ts)-(dt))/ log ((Ti-Ts)/dt) 

A = °/ 0 pi *d*L*N //uT2, heat transfer area 


51 


38 

39 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 


Q=h*A*LMTD 

printf (” the rate of heat transfer to water, is %f 
kcal/h” ,Q) 


Scilab code Exa 4.12 aniline is a tonnage oc 


//Example 4.12 

//Calculate the rise in temp 
//Given 

W=0.057 

rate of water 

W_ = W * 16.66 
di=0.0212 
Ti =32 
Tw = 80 
L = 3 

/ / Calculation 

V=(W/60) *(l/((%pi/4)*di~2)) 
//the properties of water at 
mu=7.65*10“-4 
k=0.623 

conductivity 
rho =995 
air 

cp =4.17 

heat of air 

/ / calculation 

Pr=cp*10~3*mu/k 

Re=di*V*rho/mu 

//from eq. 4.19, nusslet no. 
//from dittus boelter eq. 

Nu=0.023*Re~0.8*Pr~0.4 
f=0.0014+0.125*Re~-0.32 

//Reynold analogy 


of water . 

//m" 3/min/tube , flow 

//kg/s. water flow rate 
//m, inside diameter 
//C, inlet water temp. 

//C, wall temp . 

//m, length of pip 

//m/s, water velocity 
mean liquid temp.. 

//m" 2/s, viscosity 
//W/m C, thermal 

//kg/m"3, density of 

//kj/kg C, specific 


//Prandtl no. 

// Reynold no . 

//Prandtl no . 
//friction factor 
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St=f/2 //Stanton no. 

Nul=Re*Pr*St 

//Prandtl analogy 

Stl=(f/2)/(l+5*(Pr-l) *sqrt (f/2)) 


Nu2 = St1*Re *Pr 

//colburn analogy 


//W/nT2 C av heat 


Nu3=Re*Pr~(1/3)*(f/2) 
h = Nu3 *k/(di) 

transfer coefficient 
//Q=W_*cp *10 3* (To — Ti) =h*A*LMTD 

A = °/ 0 pi * di *L //m'2 

def f ( ’ [x] = f (To) ’ , ’ x=W_* cp*10'3* (To—Ti)—h*A* ((To—Ti) / 
log((Tw—Ti)/(Tw—To))) ’) 

To = fsolve (1 , f) 

//Revised calculation 


Tm=(Ti+To)/2 


//C, mean liquid temp. 


//the properties of water at new mean liquid temp.. 


mul=6.2*10 ~-4 
kl=0.623 

conductivity 

rho1=991 

air 

cp1=4.17 

heat of air 
//calculation 

Prl=cpl*10~3*mul/kl 

Rel=di*V*rhol/mul 

//from dittus boelter eq. 

f1=0.0014 + 0.125 * Re 1“(-0.32) 

//colburn analogy 


//m'2/s, viscosity 
//W/m C, thermal 

//kg/m'3, density of 

//kj/kg C, specific 


//Prandtl no . 

// Reynold no . 

// friction factor 


Nu4=Rel*Prl~(l/3)*(fl/2) 

hi=Nu4*kl / (di ) //W/nT2 C av heat 

transfer coefficient 

deff ( ’ [ x] = f ( To_) ’ , ’x=W_* cp*10'3*( To_—Ti)—hi *A* ((To_— 
Ti) / log ((Tw—Ti) / (Tw—To_) )) ’) 

To_=fsolve(l,f) 

printf (” Outlet temp, of water for one pass through 
the tubes is %f C”,To_) 
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Chapter 5 

free convection 


Scilab code Exa 5.1 Rate of heat loss 


//Example 5.1 

//Calculate the rate of heat 


//Given 

T1=65 //C, 

T2 =25 //C, 

h = 1.5 / / m, 

w = l / / m, 

Tf=(Tl+T2)/2 //c, 


loss . 

furnace temp, 
ambient temp . 
height of door 
width of door 
average air film temp. 


//Properties of air at Tf 


Pr = 0.695 
mu=l.85*10“-5 
beeta=l/(Tf+273) 

volumetric expension 
k=0.028 

conductivity 


// Prandtl no . 

//m A 2/s, viscosity 
//K A — 1. coefficient of 

/ /W/m C, thermal 


g = 9 • 8 


//m/s A 2, gravitational 


constant 

Grl=g*beeta*(T1-T2)*h~3/(mu~2) //Grashof no. 

Ral = Grl*Pr //Rayleigh no. 

//Nusslet no. 

Nul=(0.825+CO.387*(Ral)“(1/6))/(l+(0.492/Pr)“(9/16)) 
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~(8/27))~ 2 
hav=Nul*k/h 
coefficient 
Ad=h*w 
dt=Tl-T2 
q=hav*Ad*dt 

printf(”The rate of heat 


//average heat transfer 

//m"2 , door area 
//temp, driving force 
//W, rate of heat loss 
loss is %f W’ ,q) 


Scilab code Exa 5.2 steady state temprature 


//Example 5.2 
//Calculate the steady 
//Given 

T1 =60 
T2 =25 
h = 1 
w = 1 
q = 170 

Tf=(Tl+T2)/2 

//Properties of air at 

Pr = 0.7 

mu=l.85*10~-5 
beeta=l/(Tf+273) 

volumetric e x p e n s i o n 
k=0.028 

conductivity 
g = 9 ■ 8 

constant 


state temp, of the plate. 

//C, plate temp. 

//C, ambient temp. 

//m, width of door 
//W, rate of heat transfer 
//c , average air film temp 
Tf 

/ / Prandtl no . 

/ /nP 2 / s , viscosity 
//K" — 1. coefficient of 

/ /W/m C, thermal 

//m/s"2, gravitational 


/ / Calculation 

A = h*w 
P = 2*(h+w) 

L = A/P 


//m' 2 , plate area 

//m, perimeter of plate 

//m characteristic length 


Grl=g*beeta*(Tl-T2)*L~3/(mu"2) 
Ral=Grl*Pr 


/ / Grashof no . 
// Rayleigh no 
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//Nusslet no. 

Nul=0.54*(Ral) “(1/4) //Nusslet no. 

hav=Nul*k/L //average heat 

transfer coefficient 

Ts=q/(hav*A)+T2 

printf(”the steady state temp, of the plate is %f C” 
, Ts) 


Scilab code Exa 5.3 CALCULATE TIME REQUIRED 

//Example 5.3 

// Calculate the time required for cooling of the rod 


//Given 

d=0.0254 
1=0.4 
T1 =80 
T2 =30 
T3 = 35 
rho =7800 
rod 

cp = 0.473 

/ / Calculation 

m = y o pi/4*d~2*l*rho 
A = °/ 0 pi *d*l 
dt=Tl-T2 

difference 

h = l.32*(dt/d) “0.25 

coefficient 

i = integrate(’l/(T"(5/4)) ’ 
t=i/(3.306*A/(m*cp*10~3) ) 


//m, diameter of steel rod 
//m, length of rod 
//C, initial temp . 

//C, ambient temp. 

//c , temp, after cooling 
//kg/nU3 , density of steel 

//kj/kg C. specific heat 

//kg. mass of cylinder 
/ / m “ 2 , area of cylinder 
//c , instantaneous temp. 

//W/rrU2 C, heat transfer 

, ’T’ ,5,50) 


printfCThe required time for cooling is %f hr”,t 
/3600) 
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Scilab code Exa 5.4 Rate of heat loss 


//Example 5.4 

//Calculate the rate of heat 
per meter length of pipe. 
// given 

id=78*10~-3 
diameter 
od = 89 * 10 ~-3 
Pg = 15 

pressure 
t =2 * 10 ~-2 

preformed mineral fibre 

k=0.05 

conductivity 

Ta = 25 


loss by free convection 

/ / m, internal 

//nr, outer diameter 
//kg/cm"2, gauge 

/ / nr, thickness of 

//W/nr C. thermal 

/ / C, ambient air tenrp 


//Prandtl no 


Pr = 0.705 

// assume 

Ts =50 //C, 

1 = 1 //nr, 

Ti=200.5 //C, 

rs=od/2+t //m, 

insulation 

ri=od/2 // m , inner radius of 

insulation 

//Rate of heat transfer 
length of pipe 

Q=2*°/ 0 pi*l*k*(Ti-Ts)/(log(rs/ri)) 

//properties of air at taken at the mean 

Tf = (Ta + Ts)/2 //C 

mu = 1.76*10~-5 //nr"2/s. viscosity 

beeta=(l/(Tf+273)) //K" — 1, coefficient 

volumetric expansion 


skin tenrp. 
length 

initial tenrp . 
outer radius of 


through insulation per meter 


//W 

film tenrp. 


of 
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kl =0.027 //W/m C, thermal 

conductivity 

ds=2*rs //m, outer dia. of 

insulated pipe 

g=9.8 //m/s"2, gravitational 

constant 

Grd = g*beeta* ( Ts-Ta)*ds ~ 3/ ( mu “ 2) //Grashof no. 

Rad = Grd*Pr //Rayleigh no 


// from eq . 5.9 

//Nusslet no. 


Nu=( 0.60+(0.387*(Rad)"(1/6))/(l+(0.559/Pr)"(9/16)) 

"(8/27))"2 

hav = Nu*kl/ds //W/ in"2 C, average 

heat transfer coefficient 
Ts = (Q/(%pi*ds* 1*hav))+Ta //C, skin temp, 

//revised calculation by assuming 
Tsl=70 //C, skin temp. 

//Rate of heat transfer through insulation 
Ql=2*%pi*l*k*(Ti--Tsl)/(log(rs/ri)) 

Tf1=(Ta+Tsl)/2 //C, average aie mean 

film temp. 


mul = 1.8*10 ~-5 
beetal=(l/(Tf1+273)) 

of volumetric expansion 

kl=0.0275 


//uT2/s. viscosity 
//K" — 1, coefficient 

//W/m C, thermal 


conductivity 

Prl=0.703 //Prandtl no. 

Grdl=g*beetal*(Tsl-Ta)*ds"3/(mul"2) //Grashof 

no . 

Rad=Grdl*Prl // 

Rayleigh no. 

// from eq . 5.9 

// average heat transfer coefficient , in //W/ m A 2 C, 

havl=(0.60+(0.387*(Rad)"(l/6))/(l+(0.559/Pr)"(9/16)) 
"(8/27))"2*(kl/ds) 


Ts2 = (Q1/(%pi* ds * 1 *havl))+Ta 

//again assume skin temp. = 74 
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49Ts2=74 //C, assumed skin temp. 

50 Q3=2*%pi*l*k*(Ti-Ts2)/(log(rs/ri)) 

51 printf(”the rate of heat loss by free convection per 

meter length of pipe, is %fW’,Q3) 
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Scilab code Exa 5.5 thickness of insulation 

//Example 5.5 

//Calculate , what thickness of insulation should be 
used 

//so that the insulation skin temp, does not exceed 
65 C 


//Given 

Ts = 65 
To =30 
Tw = 460 

Tf=(Ts+To)/2 

beeta=(l/(Tf+273)) 

coefficient of volumetric 
g = 9 • 8 

constant 

mu=l.84*10“-5 
L = 10.5 
di =4 

converter 

Pr = 0.705 
k=0.0241 

conductivity 


//C, skin temp . 

//C, ambient temp. 

//C, wall temp. 

//C,mean air film temp 

//K~ — 1, 

expansion 

//m/s "2, gravitational 

//uT2/s, viscosity 
//m, height of converter 
//m, diameter of 

/ / Prandtl no . 

//kcal/h m C, thermal 


/ / Calculation 

Grl = g* beet a* (Ts-To) *L “ 3/(mu “2) //Grashof no. 

x=di/L //assume di/l=x 

y = 35/(Grl) “ (1/4) //assume 35/( Grl) “ (3/4) 

y 
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21 // printf ”x>y”” 

22 //for a verticla flat plate, from eq. 5.3 

23 Ral = Grl*Pr //Rayleigh no. 

24 // nusslet no . 

25 Nu=(0.825+(0.387*(Ral)~(l/6))/(l+(0.496/Pr)~(9/16)) 

~(8/27))~2 

26 hav=Nu*k/L //kcal/h nT2 C, average 

heat transfer coefficient 

27 //w=poly (0 ,” w” ) 

28 //Dav = (4+(4+2*w) )/2 //average 

diameter 

29 //Aav=%pi*Dav*L //average heat 

transfer area 

30 //Qi=%pi*Dav*L*0.0602* (Tw— Ts) /w //Rate of heat 

transfer through insulation 

31 //rate of heat transfer from the outer surface of 

the insulation by free convection 

32 //Qc=hav*%pi*Dav*L* (Ts—'To) 

33 //Qi=Qc 

34 deff(’[x]=f(w)’,’ x=%pi* (4+w) *L*0.0602* (Tw—Ts) /w—hav* 

%pi*(4 + 2*w) *L* (Ts— 1 To) ’) 

35 w = fsolve(0.1,f) 

36 printf (’’The required insulation thickness is %f m” , w 

) 


Scilab code Exa 5.6 rate of heat gain 
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//Example 5.6 

//Calculate the rate of heat gain by the cooler 
surface. 

//Given 


L = 1.6 
w=0.04 
b = 0.8 
T1 =22 


//m, height of enclosure 
//nr, width of enclosure 
/ / nr, breath 
//C, surface tenrp . 
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T2 =30 

Tm=(Tl+T2)/2 
Pr = 0.7 

//fpr air at 26 C 

beeta=l/(Tm+273) 

volumetric e x p e n s i o n 

mu=l.684*10“-5 
k=0.026 

alpha=2.21*10“-5 
g = 9.8 


//C, wall temp. 

//C, Mean air temp. 

/ / Prandtl no . 

//K“ — 1. coefficient of 

//m“2/s, viscosity 
//W/m C, thermal conductivity 
//m“2/s, thermal diffusity 
//m/s “2 , gravitational 


constant 

17 Raw=g*beeta*(T2-Tl)*w“3/(mu*alpha) 

no . 

18 Nuw=0.42*(Raw)“0.25*Pr“0.012*(L/w) 

no . 

19 h = Nuw*k/w 

"2 C, heat transfer coefficient 

20 q=h*(T2-Tl)*(L*b) 

rate of heat transfer 


//Rayleigh 
-0.3 // Nusslet 

//kcal/h m 
//W, the 


21 printf(”the rate of heat transfer is %fW’,q) 


Scilab code Exa 5.7 Rate of heat loss 


1 // example 5.7 

2 //Calculate the rate of heat loss by the combined 

free and forced convection . 


3 

4 

//Given 

Ts =60 

//c, 

surface temp 

5 

o 

CO 

ii 

o 

H 

//c, 

bulk temp . 

6 

d=0.06 

// m , 

diameter of pipe 

7 

1 = 1 

// m , 

length 


8 Tm=(Ts+To)/2 

9 //for air at Tin 

10 rho = 1.105 

11 cp = 0.24 
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//kg/m“3, density 
//kcal/kg C. specific 



heat 


12 mu = l . 95* 10~ -5 

13 P = 0.7 

14 kv = l .85*10~-5 

viscosity 

15 k = 0.0241 

conductivity 

16 beet a= (1 / (Tm + 273) ) 

coefficient of volumetric 

17 V = 0.3 

18 g = 9.8 


//kg/m s. viscosity 
//Prandtl no. 

//nT 2/s, kinetic 

/ / k c a 1 / f m C, thermal 

//K‘-l. 

expension 
/ / m/ s , velocity 
//m/s "2, gravitational 


constant 

19 //Calculation of nusslet no. 

20 Rad = g* beet a* (Ts-To) *d ~3* P/(kv ~ 2) //Rayleigh no. 

21 // from eq. 5.9 

22 Nufree=(0.60+(0.387*Rad~(1/6))/(l+(0.559/P)“(9/16)) 

“(8/27))~2 


23 //calculation of forced convection nusslet no. 

24 // from eq. 4.19 

25 Re=d*V/(kv) 

26 Nuforced=0. 3+(0.62*Re“(l/2)*P“(l/3)/(l+(0.4/P)“(2/3) 

)~ (1/4))*(l+(Re/(2.82*10~5))“(5/8))~(4/5) 

27 Nu=(Nuforced~3+Nufree~3)~(l/3) //nusslet no. 

for mixed convection 

28 //Nu=h*d/k 

29 h = Nu*k/d //kcal/h m' 2 C, heat 

transfer corfficient 


30 q=h* °/ 0 pi *d* 1 * ( Ts-To ) 

31 printf (” the rate of heat loss per meter length is %f 

kcal/h”,q) 
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Chapter 6 

Boiling and condensation 


Scilab code Exa 6.1 Consider nucleate pool 


1 //Example 6.1 

2 //calculate (a) the diameter of cavity on the 

boiling surface 

3 //which produce a bubble nucleus that does not 

collapse . 

4 //(b) what degree of superheat is necessary so that 

a bubble nucleus grow 

5 //in size after detachment from the cavity. 


7 Tsat=350 //K, saturated temp. 

8 Tl=Tsat+5 //K, liquid temp. 

9 //By antoine eqn. 

10 T = T1 -273 //C, 

11 pi= exp (4.22658-(1244.95/CT+217.88))) 

12 ST = 26.29-0 . 1161 * T //dyne/cm, Surface tension of 

liquid 

13 ST_ = ST*10~-3 //N/m Surface tension of liquid 

14 Lv = 33605 //kj/kgmol, molar heat of 

vaporization 


15 R = 0.08314 //m'3 bar/kgmol K, gas costant 

16 r=(2*ST_*R*Tsat~2)/ ( ( T1 - Tsat ) *pl * ( Lv * 10 ~ 3) ) 
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printf(”So a bubble nucleus that has been detached 
from a cavity will not collapse in the liquid if 
it is larger than %f micrometer \n”,r*10“6) 


//(b) 

rl=10~-6 //m 

// pil=exp (4.226 58 -(1244.95/ (Tl_ -273 + 217.88) ) ) // 

vapour pressure 

//ST1 =0.02629-1.161 *1(T - 4(T1_- 273) // 

surface tension 

deff ( ’ [x]=f (Tl) ’ , ’ x=(Tl—Tsat) 

— 2*(0.02629 —1.161*10" —4*(T1—273)) *R* Tsat "2/(rl*Lv 

* 10 " 3 ) ’) 

Tl = fsolve (0.1,f ) 

T_=(T1-273.5)-(Tsat-273) 

pr intf (” The superheat of the liquid is %f C”, round ( 
T_) ) 


Scilab code Exa 6.2 rate of boiling of water 


//Example 6.2 
//Calculate the rate of 
//Given 

d=0.35 

p=l.013 

Tl=115 

T2 = 100 

Te=Tl-T2 

//For Water 

mul=2.70*10~-4 

cp1=4.22 

rho1=958 

Lvl=2257 

vaporization 


boiling of water . 

//nr, diameter of pan 
//bar , pressure 
//C, bottom temp. 

//C, boiling temp. 

//C, excess temp. 

//Ns/m"2, viscosity 
//kj/kg C, specific heat 
/ / kg / nr63 . density 
//kj/kg, enthalpy of 
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14 s 1 =0.059 

15 Pr 1 = 1.76 

16 //For saturated steam 

17 rho2 =0.5955 

18 //For the pan 

19 Csf =0.013 

20 n = l 

21 g = 9.8 


//N/m , surface tension 
// Prandtl no . 


// constant 
/ / exponent 

//m/ s “ 2, gravitational 


constant 


22 //from eq. 6.6 //heat flux 

23 Qsl=mul*Lvl*(g*(rhol-rho2)/si)~(1/2)*(cpl*Te/(Csf* 

Lvl*(Prl)~n))~3 


//kg/m'2 s. rate of boiling 
//m~2, pan area 
//kg/s, Total rate of 


24 Rate=Qsl/Lvl 

25 Ap = °/„pi/4*d~2 

26 Trate=Rate*Ap 

boiling 

27 Trate_=Trate*3600.5 

boiling 


//kg/h. Total rate of 


28 printf (” total rate of boiling of water is %f kg/h \n 

”,Trate_) 

29 

30 //using Lienhard ’ s eq. , //critical heat flux 

31 Qmax = 0.149*Lvl*rho2*(sl*g*(rhol-rho2)/(rho2)~2) 

“(1/4) 

32 //by Mostinski eq. 

33 P c = 2 21.2 //critical pressure 

34 Pr=p/Pc //reduced pressure 

35 hb=0.00341*(Pc)~(2.3)*Te~(2.33)*Pr~(0.566) // 

boiling heat transfer coefficient 

36 hb_ = hb/1000 //kW/m“2 C boiling heat 

transfer coefficient 

37 Qs2 = hb_ * (Te) 

38 printf (”Qs2 compares reasonably well with the Qsl”) 


Scilab code Exa 6.3 formaldehyde is one of 


66 



1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 


//Example 6.3 

//Calculate the rate of boiling. 

//Given 

A = 12.5673 
B=4234.6 
pv=l.813 

Tl=200 //C, tube wall temp. 

//For methanol 

Tc=512.6 //K, critical temp. 

w = 0.556 //acentric factor 

Zra=0.29056-0.08775*w 

R=0.08314 //nL3bar/gmol K, universal gas 

constant 

Pc =80.9 //bar, critical temp. 

Mw = 32 //g , molecular wt 


/ / Calculation 

//Estimation of liquid and 
//from antoine eq. 

T=B/ (A-log (pv)) //K, 

Te=(Tl+273)-T //K, 

Tm=((Tl+273)+T)/2 //K, 


vapour properties 

boiling point 
excess temp . 
mean temp . 


//Liquid properties 

//(a) 

Tr=T/Tc //K, reduced temp, 

//from Rackett technique 


//m~3/kg mol, 
//kg/nR3 , 


Vm=R*Tc*(Zra)“(l+(l-Tr)~(2/7))/Pc 

molar volume 
rhol=Mw/Vm 

density of satorated liquid density 

//(b) 

//from Missenard technique 

T2=348 //K, given data temp. 

T3 = 373 //K, given data temp. 

Cp2 = 107.5 //j /g m °l K specific heat at T2 

Cp3 = 119.4 / / j / g m °l K specific heat at T3 

//By linear interpolation at T = 353.7 K 
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Cp=Cp2+(Cp3-Cp2)*((T-T2)/(T3-T2)) //kj/kg mol C, 

specific heat at T = 353.7 K 
Cp_=Cp*0.03125 //kj/kg C 

//(c) Surface tension at given temp . (K) 

T4 = 313 
St4 =20.96 
T5 =333 
St5 = 19.4 

//By linear interpolation at T = 353.7 K 

S=17.8 //dyne/cm, 

surface temp. 

/ / (d) liquid viscosity 
T6 =298 

MUt 6 =0.55 //cP, liquid 

viscosity at temp =298 

MU=((MUt6)“-0.2661+((T-T6)/233))“(-I/O.2661) 

//cP 

//(e)Prandtl no. a,b,c are constant 

a=0.3225 

b=-4.785*10“-4 

c = 1.168*10“-7 

kl = a+b*T + c*T“2 //W/m C, thermal 

conductivity 

Prl=Cp_*1000*MU*10“-3/kl //Prandtl no. 

//(f)heat of vaporization at 337.5 K 

Lv = l 100 //kj/kg, enthalpy 

of vaporization 

//Properties of methanol vapour at Tin 

//(a) 

Vml =R*Tm/pv //m“3/kg mol, molar 

volume 

rhov = Mw/Vml //kg/m“3, density 

of vapour 

//(b) al,bl,cl,dl are costants 

al=-7.797*10“-3 
bl=4.167*10“-5 
c1=1.214*10“-7 
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dl=-5.184*10“-11 

//thermal conductivity of vapour 

kv=al+bl*Tm+c1*Tm“2+dl*Tm“3 //W/m C 

//(c)heat capacity of vapour, a2 , b2 , c2 , d2 are 
costants 

a2 =21. 15 
b2 =7.092*10“-2 
c2 =2.589*10“-5 
d2=-2.852*10“-8 

//heat capacity of vapour, in kj/kh mol K 

Cpv=a2+b2*Tm+c2*Tm“2+d2*Tm“3 


//(d) viscosity of vapour 

T7 = 67 
MUt7 = 112 
T8 = 127 
MUt8 = 132 

//from linear inter polation at Tin 

MUv=l.364*10“-5 //kg/rn s 


// from Rohsenow ’ s eq . 

Csf=0.027 //constant 

n=1.7 //exponent value 

// from eq . 6.6 

g=9.8 //m/s~2, gravitational 

constant 

//heat flux / /kW/m“ 2 

Q=MU*10“-3*Lv*(g*(rhol-rhov)/S*10~-3)“(1/2)*(Cp_*Te 
/(Csf*Lv*(Prl)“n))“3 
// from eq . 6.11 

//from eq 6.11, critical heat flux 

Qmax = 0.131*Lv*(rhov)“(1/2)*(S*10“-3*g*(rhol-rhov)) 
“(1/4) 

/ / dimensionless radius r _ 

r=0.016 

r_ = r*(g*(rhol-rhov)/(S*10~-3)) “(1/2) 

//peak heat flux 

Qmaxl=Qmax*(0.89+2.27*exp(-3.44*sqrt(r_))) 
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// from eq . 6.12 

//heat transfer coefficient hb 

d = 0.032 //m, tube diameter 

hb=0.62*((kv"3)*rhov*(rhol-rhov)*g*(Lv*10~3+0.4*Cpv* 
Te)/(d*MUv*Te))~(1/4) 

Qb = hb*Te //kw/m"2, heat flux 

BR = Qb * 10 "-3/Lv //kg/m"2s, boilng rate 

printf (” The boilins rate is %f kg/m"2 h”,BR*3600) 


Scilab code Exa 6.4 A mixture of benzene 
//Example 6.4 

//Calculate the physical properties of the liquid. 
//Given 

Wl=200 //kg/h, rate of entering toluene 

muv=10~-5 //kg/m s, viscosity of toluene 


vapour 

mul=2.31*10"-4 
rhol=753 
rhov =3.7 

vapour 


//kg/m s, viscosity of benzene 
//kg/uT3, density of benzene 
//kg/m" 3, density of toluene 


Cpl=1968 


// j /kg C, specific heat of 


benzene 

kl=0.112 


//W/m C, thermal conductivity of 


benzene 

T1=160 
T2 = 120 
Te=Tl-T2 
Lv=3.63*10"5 
s=l.66*10"-2 
/ / Calculation 
w=0.125 
d=0.0211 


of he & hb 


//C tube wall temp. 

//C , saturated temp. 

//C, excess temp. 

// j /kg , enthalpy of vaporization 
//N/m, surface tension 


//m, mean step size 

//, internal diameter of tube 


G = W1/ (3600*°/„pi/4*(d"2) ) 


//kg/m"2 s, mass 
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flow rate 

Re 1 =G* (1 - w) *d/mul //Reynold no. 

Prl=Cpl*mul/kl //Prandtl no. 

/ / from eq . 6.2 3 

x=(w/(1-w))~(0.9)*(rhol/rhov)~ (0.5)*(muv/mul)~0 .1 

//let x=l/ s ucc ep s i b i 1 i t y 
// from eq . 6.22 


F = 2.35*(x + 0.231)~0.736 

liquid only reynold no.’ 


//factor signifies ’ 
to a two phase reynold 


no . 

// from eq . 7.21 

Re2 = 10 ~-4* Re 1 *F ~ 1.25 //Reynold no. 

// from eq . 6.18 

S=(l+0.12*Re2~1.14)~-l //boiling supression 

factor 

// from eq . 6.15 

hc=0.023*Rel~(0.8)*Prl~(0.4)*(kl/d)*F //W/nT2 C, 


forced convection boiling part 
// from eq . 6.16 

mulv = (l/rhov)-(l/rhol) //m~ 3/kg, 

kinetic viscosity of liquid vpaour 
dpsat=Te*Lv/((T2+273)*raulv) //N/m'2, change 

in saturated presssure 
//nucleate boiling part hb 

hb=1.218*10“-3*(kl“0.79 * Cpl~ 0.45*rhol~0.49*Te~0.24* 
dpsat"0.75*S/(s~0.5*mul~0.29*Lv''0.24*rhov''0.24)) 
h=hc+hb //W/uT2 C, total 


heat transfer coefficient 


//calculation of required heat transfer area 
a=5 //%, persentage 

change in rate of vaporization 
W2=Wl*a/100 //kg/h, rate of 


vaporization 

W2_=W2/3600 
Q = W2_ * Lv 
A=Q/(h*Te) 

transfer 


//kg/s 

//W, heat load 
//uR2, area of heat 
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//m, required 


46 l=A/(%pi*d) 

length of tube 

47 //from table 6.2 

48 T1 =0.393 

49 printf (” The total tube length is %f m”,Tl) 


Scilab code Exa 6.5 Saturated vapour pressure 

1 //Example 6.5 

2 //Calculate the rate of condensation of propane. 

3 //GIVEN 

4 rhol=483 //kg/m'3, density 

of liquid propane 

5 mul=9.1*10“-5 //P ,viscosity of 

liquid propane 

6 kl=0.09 //W/m K, thermal 

conductivity of liquid propane 

7 Lv = 326 //kj/kg. enthalpy 

of vaporization 

8 Cpl=2.61 //kj/kg K, specific 

heat of liquid propane 

9 T1=32 

10 T2=25 //C, surface temp. 

11 pi = 11.2 

12 rhov = 24.7 //kg/nT3, density of 

vapour 

13 g = 9.8 

14 h = 0.3 

15 // Calculation 

16 Lvl=Lv+0.68*Cpl*(Tl-T2) 

17 //h = 0 .943* (g*Lvl* 10' 3* rhoi *(rhol — rhov ) * kl ' 3 / (mul*L* ( 

T1-T2))) ' (1/4) 

18 / /Q=h*(L*1) * (T1—T2) 

19 //m=Q/ ( Lvl * 10'3) =1.867*10' -2*L' (3/4) 

20 Ref =30 
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//from the relation 4*m/mu=Re 

L=(Ref*mul/(4*l.867*10“-2))“(4/3) 

m=1.867*10~-2*L~(3/4) //rate of condensation 

for laminar flow 
/ / from eq . 6.3 2 

/ /Nul=h_/ kl * ( mul "2/(rhol*(rhol— rhov)*g))"(l/3) =Ref 
/ (1.08 * ( Ref) "(1.22) -5.2) 

Lp = h-L //length of plate over which flow is 

wavy 

A = Lp*l //m' 2 area of condensation 


h_=poly (0,”h_”) 

//Rate of condensation over total length^n( laminar ) + 
m( wavy) 

m2=m+h_*A*(Tl-T2)/(Lvl*10~3) 

Ref 1=4* m2/mul 

def f ( ’ [x] = f (hi) ’ , ’x=hl / kl * (mul "2/(rhol*( rhol —rhov ) *g 
) ) "(1/3) -(29.76 + 0.262* hi) / (1.08* (29.76 + 0.262* hi) 

" (1.22) -5.2) ’) 

hl = fsolve (1000 , f ) //W/m"2C 

m2=m+hl*A*(Tl-T2)/(Lvl*10~3) 

Ref 1=4* m2/mul 

m2=m+hl*A*(T1-T2)/(Lvl*10~3) 

printf ( r Total rate of condensation is %f kg/h”,m2 
*3600) 


Scilab code Exa 6.6 Trichloro ethylene 
//Example 6.6 

//Calculate the rate of condensation of TCE 
//(a) on a single horizontal tube 
//(b) in a condenser 
//Given 

//data fot TCE 
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7 

T1=87.4 

point 

//C, normal boiling 

8 

T2 =25 

//C, surface temp. 

9 

Lv = 320.8 

vaporization 

//kj/kg, heat of 

10 

cp = 1.105 

heat 

//kj/kg C, specific 

11 

mu=0.45*10“-3 

viscosity 

//P. liquid 

12 

k=0.1064 

conductivity 

/ /W/m C, thermal 

13 

rhol=1375 

density 

//kg/m" 3 , liquid 

14 

rhov =4.44 

vapour 

//kg/m"3, density of 

15 

Tm=(Tl+T2)/2 

//C, mean film temp . 

16 

d=0.0254 

diameter of tube 

//nr, outside 

17 

1 = 0.7 

/ / nr, length 

18 

g = 9 • 8 

gravitational constant 

//nr/s "2 , 

19 

/ / Calculation 


20 

/ / (a) from eq . 6.34 


21 

Lvl=Lv+0.68*cp*(Tl-T2) 


22 

h=0.728*(g*Lvl*10~3*rhol*(rhol- 
T2)))“(1/4) 

rhov)*k~3/(mu*d*(Tl- 

23 

A = °/„pi *d*l 

//nr"2 , area of tube 

24 

Q = h* A *(T1-T2) 

transfer 

//W, rate of heat 

25 

ra=(Q/Lvl)/1000 

condensation 

//kg/s rate of 

26 

27 

printf(”Rate of condensation is 

%f kg/h \n”,m*3600) 

28 

//(b) from eq. 6.3 5 


29 

N = 6 

vertical tire 

/ /No. of tubes in 

30 

hl=0.728*(g*Lvl*10~3*rhol*(rhol 
T1-T2)))“ (1/4) 

-rhov)*k~3/(N*mu*d*( 
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TN = 36 

TA = TN * 7 0 pi *d* 1 
Ql=hl*TA*(Tl-T2) 

transfer 


//total no. of tubes 
//uT2, total area 
//W, rate of heat 


ml=(Q1/Lvl)/1000 //kg/s rate of 

condensation 

printf(”Rate of condensation is %f kg/h \n \n”,ml 
*3600) 


//from chail ’s corelation 

h2=(l+0.2*cp*(Tl-T2)*(N-1)/(Lvl)) 
printf(”thus there will be increase in the 

calculated rate of heat transfer and in rate of 
condensation as %f percent”,18.7) 


Scilab code Exa 6.7 Saturated vapour 


//Example 6.7 

//What fraction of vapour woll condense 


//Given 

Gv = 20 

benzene 

di=0.016 
muv=8.9*(10“-6) 
Lv=391 

vaporization 

c p1 = 1.94 
Tv = 80 

benzene 
Tw = 55 
g = 9.8 

constant 
rhol=815 
rhov =2.7 
vapour 


//kg/m“2 s, mass flow rate of 

//m, tube diameter 
/ / P , viscosity 
// kj /kg . , enthalpy of 

//kj/kg C, specific heat 
//C, normal boiling point of 

//C, wall temp. 

//m/ s “ 2 , gravitational 

//kg/uT3, density of benzene 
//kg/uT3, density of benzene 
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kl=0.13 

mu=3.81*10“-4 

1=0.5 

//W/m C, thermal conductivity 
//P, viscosity of benzene 
//m, length of tube 

/ / calculation 

Rev=di*Gv/muv 

//Reynold no. of vapour 


/ / from eq . 6.3 8 

Lvl=Lv+(3/8)*cpl*(Tv-Tw) 

//heat transfer corfficient , h 

h=0.555*(g*rhol*(rhol-rhov)*kl~3*Lvl*10~3/(di*mu*(Tv 


-Tw))) “(1/4) 

Aavl = %pi * di* 1 
Q=Aavl*h*(Tv-Tw) 
m = Q /(Lvl*10~3) 

of benzene 

//m A 2 , available area 

//W, rate of heat transfer 

//kg/s, rate of condensation 


Ratei = Gv * ( °/ 0 pi/4) * di ~2 //kg/s rate of input of 

benzene vapour 

n=m/Ratei 

printf (” fr act ion of input vapour condensed is %f” , n 
* 100 ) 
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Chapter 7 

radiation heat transfer 


Scilab code Exa 7.3 the sun may be considered 


//Example 7.3 

//calculate (a)the fraction of solar radiation falls 
in visible range 

//(b) the fraction occurs on the left of visible 
range 

//(c) the fraction ooccurs on right on visible range 
//(d) wavelength and frequency of maximum spectral 
emissive power 

//(e)the maximum spectral emissive power 
//(f)the hemispherical total emissive power 
//Given 

Ts=5780 //K, surface temp. 

//CAlculation 

//(a) 


lamdal=0.4 

spectrum range 

lamda2 =0.7 

spectrum range 
El=lamdal*Ts 
E2 = lamda2 *Ts 

// from table 7.2 


//micrometer , starting visible 

/ / micrometer , ending visible 

//micrometer K, 

//micrometer K, 
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//fraction of radiation lying between 0 and lamdal 

FI=0.1229 

//fraction of radiation lying between 0 and lamda2 

F2 = 0.4889 

//the fraction of radiation falls betweem lamdal & 
lamda 2 

F3=F2-F1 

printf (”the fraction of radiation falls in visible 
range is %f \n”,F3) 

//(b) 

F4 = F1 

printf (”the fraction of radiation on the left of 
visible range is %f \n”,F4) 

//(c) 

F5=l-F2 

printf (” the fraction in right of visible range is %f 
\n”, F5 ) 

//(d) 

//from wein ’ s displacement law 

lraax=2898/Ts 

printf (’’The maximum wavelength is %f micrometer is”, 
lraax) 

c=2.998*10~8 //m/s, speed of light 

mu=c/lmax 

printf (’'The frequency is %f s*-l\n” ,mu) 

//(e) 

// from eq . 7.4 

h = 6.6256*10~-34 / /J s planck’s constant 

k = l . 3805* 10~-23 //J/K, boltzman constant 

Eblmax = (2* 0 / 0 pi*h*c~2*(lmax*10~-6) “ - 5) /((exp(h*c/ (lmax 
*10~-6*k*Ts)))-1) 

printf (’’the maximum spectral emissive power is %f W/ 
nT2\n”,Eblmax) 

//(0 

s=5.668*10~-8 //stephen costant 

Eb=s*Ts~4 

printf (’’the hemispherical total emissive power is %f 
W/nT 2 ” , Eb ) 
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Scilab code Exa 7.4 wavelength 

//Example 7.4 

//Determine the surface temp of blackbody and 
//wavelength of maximum emission. 

//Find the range of the spectrum in which the 
wavelength falls 


//Variables declaration 

Eb=4000 //W/m sq , Total emmisive power 

s = 5.669*10~-8 //Stephen boltzman constant 


/ / Calculation 

T=(Eb/s)~0.25 //k, surface temp, of black body 

ym = 2898/T //micro meter, 

//By weins law : Max. wavelength of emmision is 
inversaly proportional 

//to temprature . and constant is 2898 micrometer. 


//Result 

printf (” Surface temp, is %f C”,T) 
printf (” wavength is %f micrometer ”,ym) 
printf (” from fig 7.1 it falls in the infrared 
region of spectrum.”) 


Scilab code Exa 7.5 spectral emissivity 
//Example 7.5 

//calculate (a) total (hemispherical) emissive power 

//(b) total (hemispherical) emissivity 

//Given 
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T=1500 //K, surface temprature 

//from fig 7.7 

el=0.2 //emissivity ,when wavelength (11 ) is 0< 11 

<2 micrometer 

e2=0.6 //emissivity ,when wavelength (12 ) is 2< 12 

<6 micrometer 

e3 = 0.1 //emissivity ,when wavelength (13 ) is 6< 13 

<10 micrometer 

e4 = 0 //emissivity ,when wavelength ( 14 ) is 14 

>10 micrometer 
// from table 7.2 

FI=0.2733 

( 11 ) 


F2 =0.89-F1 

( 12 ) 

F3 = 0.9689-0.89 

(13) 

/ / Calculation 

s = 5.669*10~-8 
Eb = s *T ~4 
E=(el*Fl+e2*F2+e3*F3)*Eb 


//fraction of energy in 
//fraction of energy in 
//fraction of energy in 


wavelength 

wavelength 

wavelength 


//Stephen’s constant 
//emissive power 


printf (” total (hemispherical) emissive power is %f W 
/nC2\n”,E) 

//(b) 

e=E/(s*T~4) 

printf (” total (hemispherical) emissivity of the 
surface is %f ” ,e) 


Scilab code Exa 7.6 fraction of radiation 
//Example 7.6 

//Calculate the fraction of radiation emitted by the 
surface. 

r i = 5 II cm ,inside radius of ring 

w = 3 //cm, width 
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ro=ri+w 
L = 20 


//cm , outside radius 
/ / cm, surface distance 


//view factor along surface dAl—A2 

FI =2* integrate ( ’ 2(T 2* r / (2(T 2+r ~2) ~2 ’ , ’r ’ ,0,ri) 

//view factor along surface dAl—A2” 

F2 =2* integrate (’20*2*r/(20*2+ r *2) '2 ’ , ’ r ’ ,ri,ro) 
printf (” fraction of radiation passes through hole 
%f \n”,FI) 

printf (”fract ion of radiation intercepted by the 
ring %f ”,F2) 


Scilab code Exa 7.8 relevant view factor 
//Example 7.8 

//Consider an enclosure consisting of a hemisphere 
//of diameter d and a flat surface 
//of the same diameter. 

//Find the relevant view factor 

//Variables declaration 

F11=0 //view factor 

d=l //let it be 

printf (” view factor Fll = %f” ,F11) 

/ / Calculation 

F12 = 1-F11 //view factor 

printf ( ” view factor F22 =%f”,F12) 

A1 = ((°/ 0 pi ) *d~2) /4 //sq m, area 

A2 = ((°/„pi ) *d~2) /2 // sq m, area 

F21 = A1/A2 //from eq . 7.26 

printf (” view factor F21 =%f” , F21) 

F22=l-F21 

//Results 
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printf (” view factor=%f” , F22) 


Scilab code Exa 7.9 determine the view factors 
// Example7.9 

//Consider an enclonsnre formed by closing one end 
//of a cylinder ( diameter= D, heigh t=H) by a flat 
surface 

//and the other end by hemispherical dome. 
//Determine the view factor of all the surfaces of 
the enclousure 

//if height is twice the diameter. 

/ / 1 ,2 ,3 ,4 are given surface of enclosure in fig. 
7.21 

//Variable declaration 
s = 3 //no . of surface 

t v f = s ~ 2 //total view factor 

//using the result of example 7.8 

F11=0 
F33= 0.5 

printf (” view factor Fll=%f”,Fll) 
printf (” view factor F33=%f”,F33) 

//Calculation & Results 
R1=0.25 //R=d/2*h &h=2d 

R2 =0.25 

X = l+((1 + R2~2)/(R1~ 2)) 

F14 = (0.5)* (X-sqrt ((X “2) -4*(R2/R1)“2)) 
printf (” view factor F14=%f”,F14) 

F13=F14 

printf (” view factor F13=%f”,F13) 

F12 = l-Fll-F13 // from eq. 7.31 for surface 1 

printf (” view factor F12=%f”,F12) 
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d=l //say 

Al = (°/ 0 pi*(d~2) ) /4 
A3 = (°/ 0 pi*(d~2) ) /2 
F31=A1*F13/(A3) 


printf (” view factor F31=%f”,F31) 

// from eq. 7.31 for surface 3 

F33 =0.5 
F32=l-F31-F33 

printf (” view factor F32=%f”,F32) 

//for surface 2 

A2=2*%pi*d~2 

F21=A1*F12/A2 

printf (” view factor F21=%f”,F21) 
F23=A3*F32/A2 

printf (” view factor F23=%f”,F23) 
F22=l-F21-F23 


printf (” view factor F22=%f”,F22) 


Scilab code Exa 7.10 view factors 
//Example 7.10 

//Calculate the view factors of the surfaces. 

//Given 

ds = 0.3 / / m, diameter of shell 

rl=0.1 //m, distance from the centre 

/ / Calculation 

//by the defination of view factor 

F12 = 1 

printf (’’The view factor from surface 1 to 2 is %f\n” 
, F12 ) 

//F21 

R=ds/2 //m, radius of sphere 

r2 = sqrt (R“2-rl~2) 
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A1 = °/„pi *r2 ~ 2 //in' 2 area 

A2=2*%pi*R~2 + 2* %pi*R* sqrt (R~2-r2~2) 

//from reciprocity relation 

F21=(A1/A2)*F12 

printfCThe view factor from surface 2 to 1 is %f\n” 
, F21) 


Scilab code Exa 7.12 a carbon steel sphere 


//Example 7.12 

//calculate the time required for ball to cool down. 
//Given 

d = 0.3 / /m, diameter of steel sphere 

Ti=800 //K, initial temp, of sphere 

T2=303 //C, ambient temp. 

Tl=343 //C, final tempreture 

rho=7801 //kg/m~3, density of steel 

cp = 0.473 //kj/kg C, specific heat of steel 
/ / calculation 

R = d/ 2 / /m, radius of sphere 

Al=4*%pi*R~2 //m~2, area of sphere 

m = 4/3* °/ 0 pi *R~3* rho //m' 3 , mass of sphere 
F12 = 1 //view factor 

s = 5.669*10~-8 //stephen Boltzman ’ s constant 

//—dTl/ dt=Al* F12* s * (T"4—T2 " 4) / (m* cp ) 

I = integrate( ’(1/( Tl~4—T2"4)) ’T1’ ,343,800) 

t=1/(Al*F12*s/(m*cp*10~3)) 

printf(”The time required for the ball to cool is %f 
h”,t/3600) 


Scilab code Exa 7.13 A schedule pipe 

//Example 7.13 
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//Calculate the net rate of heat loss 
//from unit length of pipe by radiation if 
//(a) tha pipe surface is considered black 
//(b) the pipe surface has an emissivity of 0.74 


//Variables declaration 

d=0.114 //m, dia.o f pipe 

1 = 1 //m, length of pipe 

A=(Zpi)*d*1 //m sq, area 

el = l //emmisivity of black body 

F12= 1 //view factor , l:pipe surface , 2: 

room walls 


s = 5.67*10~-8 //stephen boltzman constant 

Tl= 440 //K, steam temp. 

T2=300 //K, wall temp. 

/ / Caluclation 

Q12 = A* e1 * F 12 *s*(T1~4-T2"4) //net rate of radiative 

heat loss 


//Results 

printf(”(a) Net rate of radiative heat loss Q12 —%i 
W \n”,Q12) 

// Part—b 

e2 =0.74 

Q12=A*e2*F12*s*(T1~4-T2~4) //net rate of radiative 
heat loss 

printf(”(b) Net rate of radiative heat loss Q12 =%f 
W’ ,Q12) 


Scilab code Exa 7.14 view factors and rate of loss 
//Example 7.14 

// a. Calculate i—View factors F12 and //F21 , ii — 
Calculate net rate of radiant energy gain by 
inner surface . 
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//(b) Hence calculate the rate of loss 

//of saturated liquid nitrogen at 1 atm pressure 

//stored in a double walled spherical Dewar flask. 


//Variable 

F12 = 1 
r1=0.15 
r2 =0.155 


declaration 

/ / view factor 
//m inner radius of phere 
//m , outer radius 


/ / Calculation 

A1 =4* (°/ 0 pi ) *rl ~2 //sq m inner area 
A2 =4* ( 7 0 pi ) *r2 ~ 2 //sq m, outer area 
F21=A1/A2 


h = 200 

nitrogen 
s=5.669*10~-8 
T2 =298 
T1 =77 
e1=0.06 
e2 =0.06 


//J/g , heat of vaporization of 

// boltzman constant 

//K, temp, of outer wall 
//K, Temp, of inner wall 
II e m m i s i v i t y 
II e m m i s i v i t y 


x=((l-el)/(el*Al))+(l/(Al*F12))+((l-e2)/(e2*A2)) 
Qlnet=(s*(T2~4-Tl~4))/(x) 


// Result —a—i 

printf (” a—i ) View factor F12 = %f”,F12) 
printf (” view factor F21 = %f”,F21) 

//Result— b 

printf (”(ii) The net rate of heat gain Qlnet =%f J/s 
” ,Qlnet) 
nl = Qlnet/h 

nl=nl*3600 //g/h 

printf (”(b) Rate of nitrogen loss = %f g/h”,nl) 


Scilab code Exa 7.15 Net rate of radient heat 
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/ / Example 7.15 

//Calculate the net rate of radiant heat transfer to 
the wall. 


//Given 

x=0.15 //m, length of opening on a 

furnace 

y=0.12 //m, width of opening on a 

furnace 


x 1 =6 

yi=5 
e2 = 0.8 
T1=1400 
T2 =35 
T3 = 273 

s=5.669*10~-8 


//m, width of wall 
//m, height of wall 
//emissivity of wall 
//C, furnace temp. 

//C, wall temp . 

//C, standard temp. 

//stephen boltzman ’ s constant 


//in fig. 7.29 

11=2 //nr, 11=AF 

12 = 1.5 //nr, 12=AH 

h = 3 //nr, E=dAl 

//for the dAl—A2 pair the equation is 

Fl=(l/(2*%pi)) *((12/(sqrt (12~2+h~2))) *tanh (11/( sqrt ( 
12 ~ 2 + h ~ 2)))+(11 /(sqrt (11~2+h~2)))* tanh (12/( sqrt ( 
11~2+h~2)))) 


/ / Similarly 

//for the dAl—A3 pair the equation is 

F2 = 0.1175 


//for the dAl—A4 pair the equation is 

F3 = 0.1641 

//for the dAl—A5 pair the equation is 

F4 = 0.0992 

//view factor b/w the opening (dAl)and the wall (W) 
i s 


F5=F1+F2+F3+F4 

//Calculation of radient heat exchange 

dAl=x*y 
Aw=xl*yl 
Ebl=s*(T1+T3)~4 
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Ebw=s*(T2+T3)"4 
F6=dAl*F5/Aw 

Q=dAl*F5*e2*(Ebl*(l-(l-e2)*F6)-Ebw) 

printf(”the net rate of radiant heat transfer to the 
wall is %f W’ ,Q) 


Scilab code Exa 7.16 the base of rectangular 


//Example 7.16 

//Part—a—If the side walls are perfectly insulated 
//and the surfaces are diffuse gray 
//with an emissivity 0.7 

//, Calculate the required net rate of heat supplied 
to base . 

//b— If the skin temp, of the outside of the top 
wall is 60 degree celcius 
//and heat loss frirn this surface occurs 
//to a big factory shade at 30 degree celcius 
//calculate the convective heat transfer coefficient 


//Variable declaration 

1=3 / /:m, length of wall 

w = 2 //m, width of , wall 

d=3 //m 

Rl=l/d 


Al=l*w 
A2 = A1 
e1=0.7 
e2 =0.7 
T1=673 
T2 = 523 

s = 5.669*10~-8 


//sq m, area 1: front part 
//sq m , area, 2” back part 
/ / e m m i s i v i t y 
/ / e m m i s i v i t y 

//k 

// k 

//stephen boltzman constant 


/ / Calculation 

F12 = 0.148 //view factor , from fig. 7.12 

x=(Al+A2-2*Al*F12)/(A2-(A1*(F12~2)))+((l/el)-l)+(Al/ 

88 



26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 


A2)*((l/e2)-1) 


// Results 

Qlnet=-l*Al*(s*(T2~4-Tl~4))/(x) 

printf (’’the net rate of radiant heat loss =%f kW \n” 

,Qlnet/1000) 

// (b) 

F24 = l //from fig 7.12 

T20=333 //K, outer surface temp, of surface 2 

T4=303 //K, ambient temp 

Q2rad=A2*e2*F24*s*(T20~4-T4~4) 
q=Qlnet-Q2rad 
ql=q/1000 // Kw 

h=q/(A2*(T20-T4)) 

printf (” convective heat transfer coeff. =%f W/sq mC 
”,h) 


Scilab code Exa 7.17 two parallel disks 
//Example 7.17 


//calculate 
between 
//given 

the net rate 
the disks. 

of exchange of 

radiation 

rli=0.1 

// m , 

inner 

radius of 

disk 

1 

r1o =0.2 

// m > 

outer 

radius of 

disk 

1 

r2i=0 . 12 

// m > 

inner 

radius of 

disk 

2 

r2o =0.25 

// m > 

outer 

radius of 

disk 

2 

h=0.08 
R2=r2o/h 

// m > 

distance between 

the 

disks 


Rl=rlo/h 

X=1+(1+R1~2)/R2~2 

F23_14=l/2* (X-sqrt (X~2-4*(R1/R2)"2)) 
//calculation of F23_4 

R2_=r2o/h 

Rl_=rli/h 
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X_ = 1+(1 + R1_~2)/R2_ ~ 2 

F23_4=l/2*(X_-sqrt(X_~2-4*(Rl_/R2_)~2)) //view 

factor 
//similarly 

F3_14 = 0.815 //view factor 

F34 =0.4 //view factor 

A23 = °/ 0 pi*r2o~2 //area 

A3 = °/ 0 pi *r2i ~2 
Al = '/.pi*(rlo“2-rli“2) 

// from eq . 1 

F12=A23*(F23_14-F23_4)/A1-(A3*(F3_14-F34))/A1 


/ / calculation 
//given 

T1=1000 
T2 =300 

s = 5.669*10~-8 
el=0.8 
e2 =0.7 


of the rate of radiative heat exchange 

//K, temprature of disk 1 
//K, temprature of disk 2 
//Stephen’s Boltzman constant 
// emissivity 


A2 = °/ 0 pi*(r2o~2-r2i“2) 
Fls=l-F12 


F2s=l-(A1*F12/A2) 

/ / calculation 

//let some quantities equal to 

a=(l-el)/(el*A1) 
b=l/(Al*F12) 
c=(l-e2)/(e2*A2) 
d=l/(Al*Fls) 
e=l/(A2*F2s) 


f=s*Tl~4 


g=s*T2"4 

// from eq . 7.42(a) 

/ / (f—Jl)/a = (Jl—J2) /b+Jl / d 
//(g—J2)/c = ( J2—J1)/b+Jl/e 
//solving two eqns by matrix 

A=[-0.0564,0.5036;0.4712,-0.0564] 
B= [161.847;21376.31] 

X= inv (A)*B 
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J1=X( [1] ) 

J2 = X( [2]) 

//net rate of radiation exchange 

Q12net=(J1-J2)/b 

printf (”net rate of radiation exchange b/w disk 1 
and 2 is %f W/m* 2” , Q 12net ) 


Scilab code Exa 7.18 rate of heat gain 


di=0.0254 
Ti =77 

do =52.5*10~-3 
To = 270 
1 = 1 

e1=0.05 
e2 = 0.1 
e3 = 0.02 


//Example 7.18 

//calculate the rate of heat gain by the liquid. 
//Given 

//m, inner diameter of tube 
//K, liquid temprature 
//m, pipe internal diameter 
//K, wall temprature 
//m, length of tube 
//emissivity of tube wall 
//emissivity of pipe wall 
//emissivity for inner surface of 
radiation field 

e4 = 0.03 //emissivity for outer surface of 

radiation field 

s = 5.669*10~-8 //stephen boltzman costantl 

/ / Calculation 

ds=(do+di)/2 
Ao = %pi*do * 1 
As = 7„pi *ds * 1 
Ai = °/„pi *di *1 

//View factors 

//for the long cylindrical enclosure made up of the 
outer pipe and the shield 

Fso=l //because outer surface of shield cant see 
itself 

Fos=As/Ao 


//m, diameter of radiation shield 
//m~2, outer pipe area 
//m "2, shield area 
//m "2, inner pipe area 
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Fsi=Ai/As 

//now assume 

/ / (1 — e2)/e2+ 1/Fos +Ao*(l — e4 ) / (As* e4 )=x 
// (1 — e3)/e3 +1/Fsi +(1/Fsi) *(1 — el)/el=y 

x=(l-e2)/e2+ 1/Fos +Ao*(1-e4)/(As *e4) 
y=(l-e3)/e3 +1/Fsi +(1/Fsi)*(1-el)/el 

//solving the equations for heat transfer from the 
outer pipe and inner pipe 

deff(’[x]=f(Ts)’,’ x=(Ao* (To~4—Ts~4) /x) — (Ai * (Ts~4— Ti 
~4)/x) ’) 

Ts = f solve (1,f) 

Qos=(Ao*s*(To~4-Ts~4))/x 

printf(”The net rate of heat gain of tube is %f W’ , 
Qos) 


Scilab code Exa 7.20 carbon dioxide gas 

//Example 7.20 

//Calculate the spectral extinction coefficient. 

// Given 

T = 300 
per = 91 
lam=4.2 
L = 0 . 1 

/ / calculation 
// 12/11 = f 

f=l-per/100 

transmitted 
// from eq . 7.69 

a=- log (f)/L 

printf (” the spectral extinction coefficient is %f m 
M”,a) 


//K, temprature 
//percent , adsorbed radiation 
//micrometer, wavelength radiation 
//nr, path length 

//fraction of incident radiation 
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Scilab code Exa 7.21 hot flue gas 


//Example 7.21 
//Calculate the 
//Given 

rate of heat transfer . 


Ts = 800 

//C, wall temp. 


Tg=l100 

//C. burner temprature 


C02 =8 

//percent , composition of 

C02 in flue 

gas 



M=15.2 

//percent , composition of 

moisture in 

flue gas 



a = 0.4 

//m, length of duct 


o 

ii 

,Ci 

//width of duct 


h = 15 

//W/m"2 C, heat transfer coefficient 

P = 1 

//atm pressure 


//CA1CULATI0N of Eg(Tg) 


pc = C02 / 100 * P 

//atm, partial pressure of C02 

pw=M/100*P 

//atm, partial pressure of 

moisture 



1 = 1 

//m, length of duct 


V=a*b*l 

//m*3, volume of duct 


A = 1.6*1 

//m~2 area of duct 


Le =3.6*(V/A) 

//m, mean beam length 


pc*Le 

pw*Le 

Tg_=Tg+273 
Ts_=Ts+273 
// from fig 7.3 8 

Ec =0.06 

Eg = 0.048 

// from fig 7.3 9 


//a correction 

pw/(pc+pw) 
pc*Le+pw*Le 

dE need to be calculated 


// from fig . 7.3 9 


dE=0.003 
Eg_Tg=Ec+Eg-dE 

//emissivity at temp. 

Tg 
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//Calculation of alpha 

p c * L e * T s / T g 

/ / from fig. 7.37 
Ec1=0.068 
//from fig . 7.3 8 

Ewl=0.069 

Cc=l //correction factor 

Cw=l //correction factor 

d_alpha = dE //AT 1 ATM TOTAL PRESSURE 

alpha=Cc*Ecl*(Tg_/Ts_)~0.65+Cw*Ewl*(Tg_/Ts_)~ 0.45 -dE 

//radiant heat ransfer rate 

s = 5.669*10~-8 //Stephen’s boltzman 

constant 

Qrad = A*s*(Eg_Tg*Tg_~4-alpha*Ts_~4) / /kW 

Qconv=h*A*(Tg-Ts) //kW, convective heat transfer 

rate 

Q=Qrad+Qconv 

printf(”The total rate of heat transfer from the gas 
to the wall is %f kW” ,Q/1000) 
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Chapter 8 

Heat Exchanger 


Scilab code Exa 8.1 Benzene from condenser 


//Example 8.1 
//page no. 303 
//Given 

//for Benzene 

Mb=1000 

//Kg, mass of benzene 

T1 =75 

//C initial temp, of benzene 

T2 =50 

//C final temp, of benzene 

Cpl=l.88 

//Kj/Kg C. specific heat of 

benzene 

mul=0.37 

//cP. viscosity of benzene 

rho1=860 

//kg/m" 3, density 

kl=0.154 

//W/rnK. thermal conductivity 

//for water 

Tav=35 

//C av , temp . 

Cp2 =4.187 

//specific heat 

mu2 =0.8 

//cP. viscosity 

k2 =0.623 

//W/mK. thermal conductivity 

o 

CO 

ii 

co 

H 

//C. initial temp . 

H 

II 

4^ 

o 

//C final temp . 
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/ / Calculation 

//(a) 

HD = Mb*Cpl*(Tl-T2) //Kj/li, heat duty 

WR = HD / (Cp2*(T4-T3) ) //kg/h Water rate 

printf(”the heat duty of the exchanger is %f kj/li”, 

HD) 

printf (” the water flow rate is %f kg/h”,WR) 


//(b) 

/ /tube side 
// given 

di1=21 
do 1=25.4 
t =2.2 
kw=74.5 

the wall 
di2=41 
do2 =48 


(water) calculations 

//mm, inner diameter of inner tube 
//mm, outer dia. of inner tube 
//mm/ wall thickness 
//W/m K. thermal conductivity of 

/ / mm, inner diameter of outer pipe 
//mm, outer diameter of outer pipe 


FAl = (°/„pi/4) *(dil*10~-3) ~2 //nC2 , 

FR1=WR/1000 
vl=FRl/(FA1*3600) 

velocity 

Rel=(dil*10~-3)*vl*1000/(mu2*10~-3) 
Prl=Cp2*1000*(mu2*10~-3)/k2 
//using dittus boelter eq. 
Nul=0.023*(Rel)~(0.8)*(Pr1)~(0.3) 
hl=Nul*k2/(dil*10~-3) 

transfer coefficient 


flow area 
//m/s , 

//Reynold no. 

/ / Prandtl no . 

//nusslet no . 

//W/rrh2 C, heat 


//Outer side (benzene) calculation 

FA2 = (7„pi/4) *(di2*10~-3) ~2 - ( 7„pi/4) * ( do 1 * 10 ~-3) ~2 // 

flow area 

wp = 7opi*(di2*10~-3+dol*10~-3) // 

wettwd perimeter 

dh=4* FA2/wp // 

h y d r o 1 i c diameter 

bfr=Mb/rho1 // 
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m~ 3/h benzene flow rate 

52 v2 = bf r/(FA2*3600) // 

m/ s , velocity 

53 Re2=dh*v2*rhol/ (mul *10"-3) // 

Reynold no 

54 Pr2=Cpl*10~3*(mul*10~-3)/kl // 

Prandtl no . 

55 Nu2=0.023*(Re2)“(0.8)*(Pr2)“(0.4) // 

n u s s 1 e t no . 

56 h2 = Nu2*kl / (dh) //W/nT2 

C, heat transfer coefficient 


57 printf ( ,: heat transfer coefficient based on inside 

area is %f W/nh 2 C \n”,hl) 

58 printf ( ,: heat transfer coefficient based on outside 

area is %f W/nh 2 C \n”,h2) 

59 

60 //Calculation of clean overall heat transfer 

coefficient , outside area basis 

61 // from eq . 8.28 

62 //given 

63 1 = 1 //assume , length 

64 Ao = °/„pi*dol*10"-3*l 

65 Ai = °/ 0 pi*dil*10~-3*l 

66 Am=(dol*10~-3-dil*10~-3)*%pi*l/(log(dol*10~-3/(dil 

*10~-3)) ) 

67 

68 //overall heat transfer coefficient 

69 Uo = l/ ((l/h2) + (Ao/Am)*((dol*10“-3-dil*10“-3)/( 2 * kw)) 

+(Ao/Ai)*(1/hl)) 

70 Ui = Uo*Ao/Ai 

71 

72 //Calculation of LMTD 

73 dt 1 =T1 -T4 

74 dt2 = T2-T3 

75 LMTD = (dt 1-dt2 ) / log (dt 1/dt 2 ) //log mean temp. 

difference correction factor 

76 Q = HD * 1000/3600 / /W, heat required 

77 Ao_ =Q/( Uo *LMTD) //m~@, required area 
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//m, tube length 


len=Ao_/(%pi*dol*10~(-3)) 

necessary 


//(c) 

la = 15 //m,actual length 

Aht =(%pi*dol*10~(-3)*la) 


Udo=Q/(Aht*LMTD) //W/nP2 C, overall 

heat transfer coefficient with dirt factor 
// from eq . 8.2 

Rdo = (l/Udo)-(1/Uo) //nT2 C/W 

printf (” overall heat transfer coefficient outside 
area basis is %f W/m*2 C \n”,Uo) 
printf (” overall heat transfer coefficient inside 
area basis is %f W/m*2 C \n”,Ui) 
printf (” The fouling factor is %f m'2 C/W”,Rdo) 


Scilab code Exa 8.2 design procedure 


//Example 8.2 
//Page no. 309 


//Given 

Cp = 50 
T1=135 
T2 =40 
T3 = 30 

dtl=(Tl-T2) 
dt2=(T2-T3) 
//Properties 
rho1=840 
cpl=2.093 
T = 87.5 


//tpd , plant capacity 
//C, Temp. 

//C temp . 

//C temp . 

//C hot end temp. 

//C cold end temp, 
of ethylbenzene 
//kg/nP3, density 
//kj/kg K , specific heat 

//c 

mul = exp (-6.106 + 1353/(T+ 273)+5.112*10~-3*(T + 273) 

-4.552*10 ~-6*((T + 273) ~2)) 

kl=0.2142-(3.44*10“ -4)*(T + 273)+(1.947*10~-7)*(T + 273) 
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~ 2 


kl_=kl*0.86 
// properties 

rho2 =993 
mu2=8*10~-4 
cp2 =4.175 
k2 =0.623 
k2_=k2*0.8603 


//kcal/h m K 
of water 

//kg/m A 3, density 
//kg/m s , viscosity 
//kj/kg K , specific heat 
//W/m K, thermal conductivity 
//kcal/h nT2 K 


/ / Calculation 
// ( i ) Energy balance 

Cp = Cp* 1000/24 //kg/h, plant capacity 

Cp=2083 //approx. 

HD = Cp * cpl * dt 1 //kj/h, Heat duty 

HD_ = HD *0.238837 //kcal/h 


wfr=HD/(cp2*dt2) 


//( ii ) 

raul=mul //cP, viscosity of ethylbenzene 

kl = kl //W/m K, thermal conductivity of 

ethylbenzene 


//( iii ) 

//LMTD calculation 
LMTD=(dtl-dt2)/log(dtl/dt2) 

// assume 

Udo=350 //kcal/h m~2 C, overall 

coefficient 

A = HD_ / (Udo * LMTD ) //uT2, area required 


//( iv ) 

id = 15.7 //mm, internal diameter of tube 

od = 19 //mm, outer diameter of tube 

1 = 3000 //mm, length 

0SA = °/„pi * ( od* 10"-3) * (1 * 10“-3) //uT2. outer surface 
area 

n=A/0SA //no. of tubes 

required 

f a = n* (°/,pi /4) * ( id * 10 ~-3) ~2 //nT2, flow arae 
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lv=(wfr/1000)/(3600*fa) 


//m/s, linear velocity 


//(v) 


nl=44 //total no. of tubes that can be 

accomodated in a 10 inch shell 
np=ll //no. of tubes in each pass 


//( vi ) 

bf=0.15 //m, baffel 

// ( vii ) 

//estimation of heat transfer 
//Tube side (water) 

f a 1 = (°/ 0 p i / 4) *(id*10~-3) ~ 2 * np 
vl=(wfr/1000)/(3600*fal) 
Re=(id*10~-3)*vl*rho2/mu2 

// from fig . 8.11(a) 

jh = 85 


spacing 

coefficient 

//m" 2 , flow area 
/ /m/ s , velocity 
// Reynold no . 

//colburn factor 


//jh = (hi*di) /k *(cp «nu / k) " — 1 /3 
//assume, (cp*mu/k)=x 

hi=jh*(k2_/(id*10~-3))*( C p2* 1000*mu2/k2)~(1/3) // 

kcal /h rrT 2 C 


//shell side(organic) 

c=(25.4-od)*10~-3 

adjacent tubes 

B = bf 

p=0.0254 
Ds =0.254 

shell 

/ / from eq . 8.3 2 

As=c*B*Ds/p 
Gs=Cp/As 

rate of shell fluid 
do=od/10 
shell 

//from eq. 8.31 

Dh=4*((0.5*p*100)*(0.86*p 
°/ 0 pi *do) /2) 

Dh_ = Dh * 10 ~-2 


//m, clearance b/w 2 

//m, baffel spacing 
//m, radius of 1 tube 
/ /m, inside diameter of 

//m" 2 , flow area 
//kg/m"2 h , mass flow 

//cm, outside diameter of 

100) - ((°/„pi*(do) ~2) /8) ) / (( 

//m, hydrolic diameter 
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Rel=(Dh_*Gs)/(3600*(mul*10~-3)) 

// from fig 8.11(b) 

j hi=32 

ho = jhl*(kl_/Dh_)*((6) "(1/3)) 

// from eq . 8.28 

ratio=od/id 
Rdo =0.21*10“-3 

factor 

Rdi=0.35*10~-3 


// Reynold no . 
//colb nr n factor 

// r at io=Ao/Ai 
// outside dirt 

//inside dirt 


factor 

Udo=l/((l/ho)+Rdo+(ratio)*Rdi+(ratio)*(l/hi)) 


//SECOND TRIAL 

//estimation of heat transfer 
//Tube side (water) 

npl=12 // 

fa2 = (°/„pi/4) *(id*10~-3) ~2* npl 
v2=(wfr/1000)/(3600*fa2) 
Re2=(id*10~-3)*v2*rho2/mu2 

// from fig . 8.11(a) 

j ht=83 


coefficient 


//in" 2 , flow area 
//m/s , velocity 
//Reynold no. 

//colburn factor 


//jh = (hi*di) / k*(cp *mu / k) " — 1 /3 
// assume , ( cp*mu/k)=x 

hit=jht*(k2_/(id*10~-3))*(cp2*1000*mu2/k2)~(l/3) // 

kcal /h rrL2 C 


//shell side 

c2=(25.4-od)*10~-3 
adjacent tubes 
B2 =0 . 1 
p2 = 0.0254 
Ds2 =0.254 
of shell 
//from eq. 8.3 2 
As2=c2*B2*Ds2/p2 
Gs2=Cp/As2 

rate of shell fluid 
do2=od/10 


//m, clearance b/w 2 

//m, baffel spacing 
//m, radius of 1 tube 
//m, inside diameter 

//m" 2 , flow area 
//kg/m"2 h, mass flow 

//cm, outside diameter 
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of shell 
/ / from eq . 8.3 0 

Dh2=4*( (p2*100) ~2-((°/ 0 pi*(do2) ~ 2 ) /4) ) / ((°/ 0 pi*do2) ) 

Dh2_=Dh2*10~-2 //m, hydrolic diameter 

Re2 = (Dh2_*Gs2)/(3600*(mu1 * 10~-3)) 

// from fig 8.11(b) 

jh2=48 //colburn factor 

ho2=jh2*(kl_/Dh2_)*((6)~(l/3)) 

// from eq . 8.28 

ratio=od/id //ratio=Ao/Ai 

Rdo2 =0.21 * 10~-3 //outside dirt factor 

Rdi2=0.35*10~-3 //inside dirt factor 

Udo2=l/((l/ho2)+Rdo+(ratio)*Rdi+(ratio)*(l/hit)) 


//from eq. 8.10(a) 

tauc=(T2-T3)/(T1-T3) 

R=(T1-T2)/(T2-T3) 

Ft =0.8 

Areq = HD_/(Udo2*Ft *LMTD) 
tubes =48 
lnt =4.5 

Aavl = (°/ 0 pi*od*10~-3) *tubes*lnt 
excA=((Aavl-Areq)/Areq)*100 


//Temprature ratio 
//Temprature ratio 
//LMTD correction 
// area required 
//no. of tubes 
//length of 1 tube 
// available 
//% excess 


ftor 


area 

area 


//Pressure drop calculation 
//Tube side 
//from eq. 8.3 3 

Gt = wf r / (3600*f a2) //kg/m~2 s, mass flow 

rate of tube fluid 

n2=4 //tube passes 

fit = l //dimensionless viscosity 

ratio 


g = 9.8 //gravitational constant 

f=0.0037 //friction factor 

dpt=f*Gt~2*lnt*n2/(2*g*rho2*id*10~-3*fit) //kg/ 

m~ 2, tube side pressure drop 


//eq.8.35 
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147 dpr=4*n2*v2~2*rho2/(2*g) 

tube pressure loss 

148 dpr_ = dpr *9.801 

149 tpr = dpt + dpr 

pressure drop 

150 //shell side 

151 f s =0.052 

factor for shell 

152 bf 1=0.1 

spacing 

153 Nb = lnt /bf 1-1 


//kg/nP2, return 


//N/m A 2 

//kg/nR 2 , total 


// friction 


//m, baffel 


//no. of baffles 


154 dps = fs*(Gs2/3600) ~2*Ds*(Nb+l)/(2*g*rhol*Dh2_*fit) 

//kg/nR2, shell side pressure drop 


//N/m A 2 , shell 


155 dps_ = dps *9.81 


side pressure drop 

156 printf (” Tube side Pressure drop is %f N/n/2 \n”,dpr_ 

) 

157 printf (” Shell side Pressure drop is %f N/nR2 ”,dps_) 


Scilab code Exa 8.3 The effectiveness 

1 //Example 8.3 

2 //How will the heat teansfer rate and the exit oil 

temp . 

3 //be affected if the water flow rate is increased by 


20 % 


4 

5 // Given 

6 //for hot stream 

7 Wh = 10000 


//kg/h, Rate of leaving a 


hydro 1 ic system by the oil 


8 Cph = 0.454 


//Kcal/Kg C, specific heat 


of oil 


9 Thl=85 
10 Th2 =50 


//C initial temp, of oil 
//C final temp, of oil 
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//For cold stream 

Cpc = l //Kcal/Kg C, specific heat 

of water 

Tc2=30 //C final temp, of water 

Tcl=38 //C initial temp, of water 

//from heat balance eq. 

//kg/h, Rate of leaving a hydro lie system by the 
water 

Wc=Wh*Cph*(Thl-Th2)/(Cpc*(Tcl-Tc2)) 

//For the hot stream 

Cmin = Wh*Cph //Kcal/li C. Taking hot stream 

as min. stream 
//For cold stream 

Cmax = Wc*Cpc //Kcal/li C. Taking cold 

stream as max. stream 
Cr = Cmin/Cmax //Capacity ratio 

n=(Thl-Th2)/(Thl-Tc2) //effectiveness factor 

//From eq. 8.5 7 
//No. of transfer units 

NTU = -(1 + (Cr)“2)“-(l/2)*log(((2/n)-(l+Cr)-(l+(Cr)“2) 
“(1/2))/((2/n) -(l+Cr)+(l+(Cr) ~2)“(1/2) )) 

Ud = 400 //kcal/li m“2C , overall 

dirty heat transfer coefficient 
// from eq . 8.53 

A= (NTU* Cmin)/Ud //Area required 

//if the water rate is increased by 20 %, 

a = 20 

Wc_=Wc+(Wc*(a/100)) 

Cmax_ = Wc_ *Cpc 
Cr_=Cmin/Cmax_ 

//From eq . 8.56 

n_=2*((l+Cr_)+(l+(Cr_)“2)“(1/2)* (1+exp (-(l+(Cr_)“2) 

“(1/2)*NTU)) /(1-exp (-(l+(Cr_)“2)~(1/2)*NTU))) 

~(-l) 

Th2_=Thl -(n_* (Thl-Tc2)) 

ql=Wh*Cph*(Thl-Th2) //keal/h previous rate of heat 
transfer 
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40 q2=Wh*Cph*(Thl-Th2_) //kcal/h new rate of heat 
transfer 

4f //increase in rate of heat transfer 

42 dq=(q2-ql)/ql 

43 printf(’the heat teansfer rate will be affected by 

%f percent ”,dq*100 ) 


Scilab code Exa 8.4 Thermal design 

1 //Example 8.4 

2 //calculate the time required to heat the charge. 

3 


4 

//given 


5 

p=0.0795 

/ /m. pitch of the coil 

6 

dl=0.0525 

/ / m, coil diameter 

7 

h=l.464 

section 

//m, height of the limpet ted 

8 

d2=l.5 

//m, diameter of batch 


polymerization reactor 

9 

d3 = 0.5 

//m, diameter of agitator 

10 

rpra=150 

//speed of agitator 

11 

rho =850 

//kg/m3, density of monomer 

12 

rho1=900 

//kg/m3, density of fluid 

13 

mu=0.7*10~-3 

//poise, viscosity of monomer 

14 

mul=4*10 ~-3 

//poise, viscosity of fluid 

15 

cp = 0.45 

monomer 

//kcal/kg C, specific heat of 

16 

cpl=0.5 

fluid 

//kcal/kg C, specific heat of 

17 

k=0.15 

//kcal/h mC, thermal 


conductivity 

of monomer 

18 

kl=0.28 

//kcal/h mC, thermal 


conductivity 

of fluid 

19 

Rdi=0.0002 

vessel 

//h m2 C/kcal , fouling factor 
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Rdc =0.0002 
coil 
Tci=120 
liquid 
Tvi=25 
liquid 
Tvf=80 
liquid 


//h m2 C/kcal, fouling factor for 
//C, initial temp, of coil 
/ / C, initial temp, of vessel 
/ / C, final temp, of vessel 


/ / calculation 

a= 0 / 0 pi*d2*h //outside area of the vessel 

x = 60 //%. added of the unwetted 

area to the wetted area 
ao = ((dl + (x/100) *(p-dl))/p)*a //nr"2, effective 
outside heat transfer area of vessel 
ai=6.9 //m "2,in side heat 

transfer area of vessel 


//same as outside 
area , if 
thickness is very 
small 

//vessel side heat transfer coefficient 

Re=(d3~2*(rpm/60)*rho)/mu //reynold no. 

Pr=((cp*3600)*(mu))/k 

/ / from eq . 8.6 6 

y = 1 / / x=mu / muw= 1 

Nu = 0.74*(Re~(0.67))*(Pr~(0.33)) * (y ~ (0 .14)) // 

Nusslet no 

hi = Nu *(k/d2) // 

heat transfer coefficient 


//coil side heat transfer coefficient 

v = 1.5 //m/s, linear velocity of fluid 

f a= ((°/„pi/4) *dl ~2) //m2, flow area of coil 

fr = v*fa*3600 //m3/h , flow rate of the 

fluid 

Wc = fr*rho //kg/h , flow rate 

dh=(4*(°/„pi/8)*dl~2)/(dl + (°/ 0 pi/2)*dl) //m, hydrolic 
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diameter of limpet coil 

Rel=v*rhol*dh/mul //coil 

reynold no. 

Pr1=cpl*mul*3600/kl //prandtl 

no. of the coil fluid 

/ / from eq . 8.6 8 

d4=0.0321 //m, inside 

diameter of the tube 

Nul=0.021*(Rel-(0.85)*Prl~(0.4)*(d4/d2)-(0.1)*y 
~0.14) 

hc = Nul * (kl/dh) //coil side 

coefficient 

U = 1/((1 / hi ) + ( ai / (he * ao) )+Rdi+ Rdc ) //overall heat 

transfer corfficient 

//from eq. 8.6 3 

beeta=exp(U*ai/ (Wc*cpl)) 

Wv = 2200 //kg, mass of 

fluid vessel 

t=(beeta/(beeta-l))*((Wv*cp)/(Wc*cpl))*log((Tci-Tvi) 
/(Tci-Tvf)) 

printf(”the time required to heat the charge %f min” 
,t *60) 
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Chapter 9 

Evaporetion and Evaporators 


Scilab code Exa 9.1 single effect evaporator calculation 


//Example 9.1 
// page no.391 

//calculate the rate at which heat must 
//be supplied if evapiration occurs at 
// (i ) 1 atm pressure 
//a vaccum of 650 mm Hg 
//given data 
ro=1020 
sf=4.1 
sp=3.9 
ci =5 

cw=100-ci 
cf =40 
rate = 100 
solution 
ft =25 

/ / calculation 
/ / materiel balance 

Wf = rate*ro //Kg. feed entering 

Ms = ro*ci //Kg mass of solute 

Mw = ro*cw //kg, mass of water 


// kg/m" 3, density of feed 
//kj/kg C, specific heat of the 
//kj/kg C, specific heat of the 
//initial concentration 
// cone . of water 
//final cone. 

//m~3/day, rate of cone, of aq 
// C, feed temp . 


feed 

product 
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fc = cw/ci //kg,feed concentration 

pc = (100 - cf ) / cf // kg , product concentration 

wlwp = Ms*pc //Kg, water leaving with the product 

Ws = Mw-wlwp //kg, water evaporated 

Wp=wlwp+Ms // kg, product 

//energy balance 

rt=0 //C reference temp. 

ef = sf*(ft-rt) //kj/kg,enthlpy of the feed 

//case i 

Tp = 100 //temp, of the product (because the 

solute has a ’high molecular wt ’ the boiling pt 
elevation is neglected) 

ip = sp * (Tp-rt) //kj/kg, enthalpy of the product 

iv = 2680 //kj/kg, enthalpy of the vapour 

generated at 100 C and 1 atm pr . from the steam 

table 

//refer to fig. 9.23 

//from energy balance eq . (Wf* i f+qs=Wv* iv+Wp* ip ) 

qs=Ws*iv+Wp*ip-Wf-ef //Ww=Ws 

printf(”The rate at which heat must be supplied at 1 
atm pressure is %f kj / day\n”,qs) 

//case i i 

// 650 mm Hg vaccum = 110 mmHg pressure 
bp = 53.5 //C, boiling point of water 

ip2 = sp * (bp-rt) //kj/kg, enthalpy of the product 

es=2604 //kj/kg, enthalpy of the saturated 

steam (from steam table) 

//from energy balnce eq. 
qs2=Wp*ip+Ws*es-Wf-ef 

printf (” The rate at which heat must be supplied at a 
pressure of 600 mm Hg is %f kj/day ”,qs2) 


Scilab code Exa 9.2 SINGLE EFFECT EVAPORATOR CALCULATION 
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//Example 9.2 
//Page no. 393 

//calculae the steam requirement and the no. of 
tubes 

//if the height of the calandria is 1.5 m. 


// given 

ci = 10 
cf =40 
Wf=2000 
ft =30 
rp = 0.33 
bt1=75 
sst = 115 
1 = 1.5 
sh=0.946 
lh = 556.5 
bt2 =345 
h=2150 

transfer 


//%,initial concentration 
//%, final cone 
//kg/h, feed rate 
//C feed temp. 

//kg/cm"2, reduced pressure 
//C, boiling point temp. 

//C, saturated steam temp. 

/ / m, height of calandria 
//keal/kg C, specific heat of liquir 
//keal/kg latent heat of steam 
//K, boiling point of water 
//keal/h m"2 C, overall heat 
coefficient 


si=2000*(ci/100) //kg/h, solids in 

wi = 1800 //kg/h,wate in 

Wp=si/(cf/100) //kg/h, product out 

Wv = Wf-Wp //evaporation rate 

ef = sh*(ft-bt1) 


ip = 0 

lamda.s =529.5 //keal/kg, lamda_s=is — i 1 
bpe = (273 +bt 1)-345 //boiling point elevation 
//from eergy balance eq. 

Ws=(Wp*ip+Wv*lh-Wf*ef)/lamda_s 


q=Ws *lamda_s 
A=q/(h*(sst-btl)) 
di=0.0221 
At=%pi*l*di 
N=A/At 

printf(”The steam 


//kcal/h,rate of heat transfer 
/ / m ~ 2 

//m, inside diameter 
//m" 2, area of a single tube 
//no. of tubes 
required is %f kg/h\n”,Ws) 


printf (”No. of tube are %f”,N) 
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Scilab code Exa 9.3 SINGLE EFFECT EVAPORATION 


// Example9.3 
//calculate 

// i) the steam lr . to be used in the 
// i i ) heat transfer rate required 
/ / i i i ) the steam requirement. 

// given data 


calandria 


Wf=2000 
ci =8 
cf =40 
ft =30 
vp = 660 
ssp=8 

/ / calculation 

sr = Wf *(ci/100) 
Wp=sr/(cf/100) 
ap=760-vp 

evaporator 


//kg/h, feed rate 
//% initial cone. 

//% final cone . 

//C, feed temp. 

//mm Hg, vaccum pressure 

// bar absolute , saturated steam pr . 

//kg/h, solid rate 
//kg/h,concentrated product rate 
//mm Hg, absolute pressure in the 


bt =325 
l_s =2380 
R=8.303 
w=40 
M = 18 
W = 60 
ra =2000 

dtb=(R*bt~2*w*M)/(l_s*W*m) 


//K, boiling temp, of water 
//kj/kg, latent heat 
//gas constant 
//g , mass of solute 
//g, molecular wt of solvent 
//g , mass of the solvent 
//g, molecular wt of solute 

/ / C, boiling 


point 


elevation 

bp = bt + dtb //k, boiling point of 40% solution 

dt=70 //C, from given data flux becomes 

maximum at a temp . drop =70 C 
st = bp + dt //K, saturation temp, of steam in the 

steam chest 

Sp = 2.15 // bar, from steam table , saturation lr 
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. of steam at this temp. 

sh = 4.2 //kj/kg C, specific heat of product 

rt=0 //C reference teml . 

ef = sh * (f t-rt) // kj/kg, enthalpy of the feed 

ip = sh * (54 - rt) //kj/kg, enthalpy of the product 
iv=2607 //kj/kg, enthalpy of vapour produced 

// from eq 9.6 

Wv = 1600 //enthalpy of evaporation 

q=Wp*ip + Wv*iv-Wf *ef //kj/h, heat transfe rate 

required 

hvp = 2188 //kj/kg, heat of vaporization of 

saturated steam at 397 K 
rs = q/hvp //kg/h, rate of steam supply 

printf(”The steam pressure to be used in the 
calandria is %f bar ( abs )\n” , Sp) ; 
printf(”The heat transfer rate required is %f Kj/h\n 

”,q); 

printf(”Rate of steam supply is %f kg/h” ,rs); 


Scilab code Exa 9.4 MULTIPLE EFFECT EVAPORATION 


//Example 9.4 

//calculate the evaporator areas and the steam 
economy . 

//given 


Wf=6000 
ci =2 
cf =35 
ft =50 
ssp=2 

sep=0.0139 

effect 


//kg/h, feed rate 

//%, initial concentration 

//%, final cone . 

//C, feed temp . 

//bar abs, saturated steaarn pr . 

//bar abs, maintained temp, in second 


hl=2000 //W/uL2 K, overall heat transfer 

coeffcient in 1st effect 
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h2 = 1500 //W/nT2 K, overall heat transfer 

coefficient in 2nd effect 
cp=4.1 //kj/kg k,specific heat 

/ / calculation 

si = Wf * ( ci/100) //kg/li, solid in 
wi=5880 //kg/h , water in 

Wp=si/(cf/100) //kg/h product out 

wo = Wp* (1-cf/100) //kg/li, water out with the product 
ter=wi-wo //kg/h, total evaporation rate 


//boiling temp, in the first effect 

Tl = 120 //C, Temprature 

l_sl=2200 //kj/kg, latent heat 

T2 = 12 //C, boiling point in second effect 

l_s2=2470 // kj/kg in second effect 

tatd=Tl-T2 // C,tatd=dtl+dt2 =T1—T2 , total 
available temp, drop 
// from eq . 9.20 

//hl*dtl=h2*dt2 

//solving above two equations by matrix 

A=[1,1;2000,-1500] 

C= [108;0] 


X=inv (A)*C 


dt1=X( [ 1]) 
dt2 = X ( [2] ) 

tl=Tl-dtl //temp, of steam leaving the first 

effect 

t2 = T2-dt2 //temp, of steam leaving second 

effect 

//energy balance over the 1st effect , from eq.9.14 

rt1=t1 

ef=cp*(ft-tl) //kj/kg,enthalpy of feed 

i 1 =0 

lam_sl=2330 //kj/kg 

isl=lam_sl 
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/ /Wf* e f■+Ws* 1 _ s = (Wf-Wsl) * i 1 +Wsl * i s 1 
//substituting we get, 

/ /Wsl = 0.9442*Ws— 253.4 .( 1 ) 

//energy balance over second effect 
// from eq 9.15 

// (Wf—Wsl) * i 1+Wsl* lam_s 1 =(Wf—Wsl—Ws2) * i2+Ws2* is 2 


rt2=t2 
lam_s2=2470 
is2=lam_s2 
i2 = 0 

// substituting we get 

//Ws2 = 0.8404*Wsl + 617.5.(2) 

//ter ,Wsl+Ws2 = 5657.(3) 


//solving by matrix method 

A= [0.9442 , -1,0;0,0.8404 , -1 ;0,1 , 1] 

B = [253.4;-617.5;5657] 

X= inv (A)*B 
Ws = X ( [1]) 

Ws 1 =X ( [2] ) 

Ws2 = X( [3]) 

//evaporator area 

A1=Ws*l_s1/(hi*dt1) //for 1st effect 

A2 = Ws1 *lam_s1/(h2*dt2) //for second effect 

//revised calculation 
/ /taking 

dt1_ =48 
dt2_ =60 
T1_ = T1-dt1_ 

T2_ = T2-dt 2_ 
lsl_=2335 
ls2_ =247 0 

// energy balance over first effect gives 

/ /Wsl = 0.9422Ws— 231.8.(4) 

//energy balance over second effect gives 

/ /Ws2 = 0.8457Wsl + 579.5 .(5) 

//solving eq 3,4,5 
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P= [0.9422 , -1 ,0; 0 ,0.8457 , -1 ;0,1 , 1] 

Q= [231.8;-579.5;5657] 

Y=inv (P)*Q 
Ws_ = Y ( [1] ) 

Ws1_ =Y( [2]) 

Ws2_ =Y ( [3]) 

/ / eveporator area for 1st & 2nd effect in rrL 2 

Al_=Ws_*l_sl/(hl*dtl_) 

A2_=Wsl_*lsl_/(h2*dt2_) 

EA=(A1_+A2_)/2 
SE=(Wsl_+Ws2_)/Ws_ 

printf(”The evaporator area is %f square metre \n” 

,EA) ; 

printf (” Steam economy is %f”,SE); 


Scilab code Exa 9.5 MULTIPLE EFFECT EVAPORATION 

//Example 9.5 

//Determine the maximum no. of effects to be used, 
//given 

ssp = 3.32 //bar abs , saturated steam pr . 

rp = 0.195 // bar abs, residual pr. in the 

condenser 

tl=41 //K, sun of temp, losses because of BPE 

mt=8 //k, minimum available temp, driving 

force 

/ / calculation 

sst=410 //K, saturated steam temp. 

st=333 //K, corresponding saturation temp. 

when pressure in the last effect is 0.195 bar 
ttd = sst-st //K, total temp, difference 

atd = ttd-tl // K, available temp, drop across the 

unit 

n = atd/mt //maximum no. of effect 
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printf (’’Maximum no. of effects are %f”,n); 


Scilab code Exa 9.6 MULTIPLE EFFECT EVAPORATION 


//Example 9.6 

//Calculate the heat transfer area required 
//(assuming equal area for the three effects) 
//Rate of steam consumption, Steam economy 


// given 

fc = 9.5 //%, feed concentration 

pc=50 //%, product cone. 

ft=40 // C, feed temp. 

er = 2000 //kg NaOH/h, evaporation rate 

vp = 714 //mm Hg, vaccum pr . in last effect 

//heat transfer coefficients , W/rrU2 C 
hl=6000 //for first effect 

h2=3500 //for second effect 

h3 = 2500 //for third effect 


/ / calculatii n 

Wf=er/(fc/100) //kg/h, 2 tons NaOH per hour, feed 

rate 

Wp=er/(pc/100) //kg/h, product rate 

ter=Wf-Wp //kg/h, total evaporation rate 
/ / steam 


p = 3.3 //bar , assumed saturated 

//from steam table 


//C, temp . 

//kj/kg, latent heat 
//mm Hg, pressure in the last effect 
//C, boiling point of water 


Ts = 137 
l_s =2153 
pl=760-vp 
bp = 37 

//refer to fig. 9.24 

attd=Ts-bp //C,apparent total temp, drop 
//let assume the following evaporation rate for 
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three effects in kg/h 

evl=5600 
ev2 =5680 
ev3 =5773 

//cone, in three effects 

cl=er/(Wf-evl) 
c2=er/(Wf-evl-ev2) 
c3=0.5 //given 

//boiling point elevations in three effects in C 

bpe1=3.5 
bpe2 =8 
bpe3 =39 

attda=attd-(bpe1+bpe2+bpe3 ) //actual total temp, 
drop available 

//temp, drop in three effects 
//from eq. 9.2 3 

dtl=attda*((l/hl)/((l/hl)+(l/h2)+(l/h3))) 
dt2=attda*((l/h2)/((1/hl)+(l/h2)+(l/h3))) 
dt3=attda*((l/h3)/((1/hl)+(l/h2)+(l/h3))) 

/ / from table 9.4 

//enthalpy of solution in three effects in kj/kg 

i1=486 
i2 =385 
13=460 

//enthalpy of vapour generated for three effects 
in kj/kg 
is 1=2729 
is2 =2691 
is3 =2646 

//Enthalpy of condensate over effect 1,2,3 in kj / 
kg 

il 1 =0 

il2=519 

il3=418 

//Enthalpy balance over effect 1 
ef = 145 //kj/kg,enthalpy of feed 

//from energy balance eq. 
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/ /Wsl = 0.96Ws— 3200 .( 1 ) 

//enthalpy b a lane over effect 2 

//Ws2 = 0.9146Wsl + 922.(2) 

//enthalpy b a lane over effet 3 

//Ws3 = 1.073Ws2 + 0.0343Wsl - 722.(3) 

/ / t e r =Wsl+Ws2+Ws3 = 17053.(4) 

//Solving above four eqns by matrix 

A 

= [0.96 , -1,0,0;0,0.9146 , -1 , 0 ; 0,0.0343,1.073,-1;0,1 , 1 , 1] 


B =[3200;-922;722;17053] 

X= inv (A)*B 
Ws = X ( [1] ) 

Ws 1 =X ( [2] ) 

Ws2 = X( [3]) 

Ws3 = X ( [4] ) 

//calculation of heat transfer areas iver effect 
1, 2 ,3 

Al=Ws*l_s*10~3/(hi* dt1*3600) 

A2 = Wsl*(isl-il2)* 10~3/(h2*dt2*3600) 

A3=Ws2*(is2-il3)*10~3/(h3*dt3*3600) 

// Revised dt 

avar=(Al+A2+A3)/3 
dt1_ = (A1/avar)*dt1 
dt2_ = (A2/avar)*dt2 
dt3_=attda-dtl_-dt2_ 

/ / from table 9.5 

//enthalpy of vapour generated over effect 1,2,3 
in kj/kg 
isl_=2720 
is2_ =2685 
is3_=2646 

//enthalpy of so In on 1,2,3 in kj/kg 

i1_ =47 0 
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i2_ =380 
i3_ =460 

//enthalpy of condensate over effect 1 ,2,3 in kj / 

kg 

ill_=0 
il2_=513 
il3_=412 

//enthalpy balance ove effect 1,2,3 gives 

Ws_=8854 
Wsl_=5432 
Ws2_ =5812 
Ws3_=5809 

//revised heat transfer areas for effect 1 ,2,3 in 

m" 2 

Al_=Ws_*l_s*1000/(hl*dtl_*3600) 
A2_=Wsl_*(isl_-il2_)*10~3/(h2*dt2_*3600) 
A3_=Ws2_*(is2_-il3_)*10~3/(h3*22.5*3600) 
avar_=(Al_+A2_+A3_)/3 
SE=ter/Ws_ 

printf(”The areas are now reasonably close \n”) 
printf (” Steam Rate is % f Kg/h \n”,Ws_) 
printf (” Steam economy is %f”,SE) 


Scilab code Exa 9.7 MULTIPLE EFFECT EVAPORATION 

// Exalple 9.7 

//Calculate the increase in evaporation capacity 
attainable 

//also the % change in cost of concentrating a ton 
of feed . 

//Given 

Wf =3000 //kg/li, feed 

fc=8 //%, feed concentration 

pc=40 //% product concentration 
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si=Wf*(fc/100) 
pr=si/(40/100) 


//kg , so 1 id in 
//g/h, product rate 


ft =60 
er=Wf-pr 
cost=120000 
pi=4.5 
sept =700 
cp = 0.764 


//C, feed temp . 

//kg/h, evaporation rate 
//total cost per year 
//bar, low pressure steam 

//per ton. cost of steam 
// keal/kg, specific heat 


//from table 9.6 


eep = l 

oop=400000 
oop_=600000 
wd=300 
wh = wd *24 


//atm existing evaporator pressure 
// peryear ,other operatingcost 
//per yr , for proposed condition 
//days per year . working days 
//working hr 


//EXISTING OPERATING CONDITION 
rt=0 //C,reference temp, 

ef = eep * (f t-rt) //keal/kg, enthalpy of feed 

pt=100 //C,product temp. 

i 1 = cp * (pt-rt) //keal/kg, enthalpy of soln 

isl=639 //kcal/kg,enthalpy of vapour 

generated at 1 atm (from steam table) 
l_s=496 //kcal/kg , latent heat of steam at 4.5 

bar 

T=425 //K 

//heat balance 

Ws = ( er * i s 1 +pr * i 1 -Wf * ef ) / l_s //kg/h, steam 

required 

q=Ws*l_s //ton/ hr , heat supplied 

x=q/(T-(pt+273)) //x=Ud*A 
//hourly cost 

sc = Ws /1000* ( sept) // /peril, steam cost 

1c = 100 //per h,labour cost 

oc=oop/(wh) // per h,othe cost 

tc=sc+lc+oc //total cost 

C=tc/(Wf/1000) // per ton,cost per ton of feed 
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//PROPOSED OPERATING CONDITION 

bpl=320 //K, boiling point of liquid 

dt=T-bpl 


q_=x*dt //kcal/h,rate of heat supply 

sr = q_/l_s //steam rate ton per hr 

pt_=47 //C,product temp . 

ep = cp* (pt_-rt) //kcal/kg. enthalpy of product 

ev = 618 //kcal/kg, enthalpy of vapour 


generated 
//heat balance 

/ / 24 Wf_ —582Wsl_ = 2825000 .( 1 ) 

/ / material balance 

// 4Wf_—5Wsl_=0 .(2) 


//solving by matrix method 


a=[24,-582;4,-5] 
b=[-2825000;0] 
x_ = inv (a)*b 


Wf _ = x_( [1]) 

Ws1_ = x_ ( [2]) 
ic=(Wf_-Wf) /Wf 

printf(”The increase in evaporation capacity ic %f 

percentage \n”,ic*100) 

sr_ = Ws 1_/1000 //ton per hr ,steam rate 

//hourly cost 

sc_ = Ws 1_ * sept //steam cost 

lc_=200 //labour cost rs.200/ h 

oc_=oop_/wh // other cost 

tc_=sc_/1000+lc_+oc_ 

C_=tc_/(Wf_/1000) //cost per ton of feed 

ps=(C-C_)/C 

printf (” The percentage change in the cost of 

concentrating a ton of feed is %f percentage”, 

ps * 100) 


Scilab code Exa 9.8 Mechanical vapour compression 
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//Example 9.8 

//make a mechanical vapour recompression calculation 


// given 

q = 2200 //kj/kg heat of condensation of steam 

//from example 9.1 

Qr=2.337*10~8 //kj/day rate of heat supply 


/ / calculation 

Rate=Qr/q 
Rate_=l,062*10~5 
E=2800 


//kg/day steam supply rate 

//approximate value 

//kj/kg enthalpy of compressed 


vapour 

T = 175.7 
Ts = 121 
El=2700 


//C, temprature 

//C Saturation temprature 

//enthalpy at saturation 


temprature 

ql=T-Ts 
T1=100 
E2 =419 

x=(E-El)/(E1-E2) 


//Superheat of vapour 
//C hot water temprature 
//Enthalpy at hot water temp, 
//water supplied per kg of 


superheated steam 

S = 1.044 //steam obtained after 

desuperheating 


R1=8.925*10~4 //kg/day rate of vapour 


generation 

R2=S*R1 


//Rate of recompressed sat. 


steam 


R2_=9.318*10~4 //approximate value 

SR=Rate_-R2_ 

printf(”Make up steam required is %f kg/day” , round ( 


SR) ) 
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Chapter 10 

UNSTEADY STATE AND 
MULTIDIMENSIONAL HEAT 
CONDUCTION 


Scilab code Exa 10.8 MUMERICAL CALCULATION OF UNSTEADY 
STATE HEAT CONDUCTION 


//Example no. 10.8 
//Page no. 444 

//Calculate the bottom surface , mid plane ,t°P 
surface temperatures 
//of the slab after 4 hours 
//given 

1=0.05 


margarine slab 

ro = 990 

margarine slab 

cp = 0.55 

heat of slab 

k=0.143 

conductivity of slab 


/ /m,t hic k n e s s of 
//Kg/m" 3, density of 
//Kcal/kg C, ddpecific 
//kcal/h mC, thermal 


Ti =4 
To =25 


/ / C, initial temp 
//C, ambient temp. 


123 


12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

1 

2 

3 

4 

5 

6 

7 

8 

9 


t =4 
h = 8 

/ / calculation 

Fo=k*t/(ro*cp*l~2) 
Bi=h*l/k 

// from fig. 10.6 a 

Tcbar =0.7 

Tc=To+Tcbar*(Ti-To) 

// from fig 10.6 b 
// (T-To) /(Tc-To) =0.382 

T = 0.382*(Tc-To)+To 


// hours , time 
//kcal/h m A 2 C 

// , fourier no. 

//Biot no. 

//Tcbar = (Tc—To) /(Ti-To) 
//C, centre temp. 

//c , top surface temp. 


//again from fig. 10.6 b 

Tm = 0.842* (Tc-To)+To //, mid plane temp. 

printf(”The bottom surface temperature of given slab 
is %f C”,Tc); 

printfC'The top surface temperature of given slab is 
%f CM); 

printfC’The mid plane temperature of given slab is 
%fC”,Tm); 


Scilab code Exa 10.9 NUMERIC CALCULATION OF UNSTEADY STATE 
HEAT CONDUCTION 

// ExamplelO . 9 
//Page no. 449 

//calculate : (i) time required for the cantre — 

line temp. 

//to drop down to 200 C 

//(ii)the temp, at half radius at that moment 
// ( i i i ) the amount of heat that has been transfered 
to the liquid 

// by that time per metre length of the shaft 
//given data 

Ti=870 //C, initial temp 
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//C, ambient 
/ / C, centre line 


To =30 

temp . 

Tc =200 
temp . 

h = 2000 //W/nT2 C, 

surface heat transfer coefficient 


a=0.05 

cylinder 

k = 20 

conductivity 

ro =7800 
cp = 0.46 * 10 ~ 3 

specific heat 

/ / calculation 

// i 

Bi=h*a/k 
alpha=k/(ro*cp) 

diffusivity 

Tcbar=(Tc-To)/(Ti-To) 

centre line temp. 

// from fig 10.7 a 
fo=0.51 

a 1 p h a * t / a ~ 2 

t=fo*a~2/alpha 

// ii 

//at the half radius, r/a 

T = To +0.77*(Tc-To) 


//m, radius of 

/ /W/ m C, thermal 

//kg/m*3, density 
//j/kg C, 

//Biot no . 

/ / m' 2 / C, thermal 

// dimensionless 

//fourier no. fo = 
/ / s , time 

0.5 & Bi=5 

// from fig. 10.7 b 


// iii 

x = Bi~ 2* fo 

//for x =12.75 & Bi=5.0. fig.10.9 b gives 
//q/ qi =0.83 

qi= °/ 0 pi *a “ 2* (1) *ro * cp * (Ti-To ) / / k j , initial amount 

of heat energy 

//present in 1 m 
length of shaft 
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37 q = 0.83*qi //j , amount of heat 

transferee! 

38 printf(”(i) time required for the cantre — line temp. 

to drop down to 200 C is %f s” ,t) ; 

39 printf (” ( i i ) the temp, at half radius at that moment 

is %f C ”,T); 

40 printf (”( i i i ) the amount of heat that has been 

transfered to the liquid is %f Kj”,q*10~(-3)); 
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Chapter 11 

Boundary layer heat transfer 


Scilab code Exa 11.1 water at 25 degree cclcius 

//Example 11.1 
//page no. 478 

//a—Calculate Boundary layer thickness at x = 0.5 m 
//b— Calculate local drag coeff at x = 0.5 m 
//c—Force req to hold the plate in position 
//d— shear stress at a plane , distant t/2 from the 
surface at x = 0.5 m 
//Variable declaration 

v =1 / /ill/ s 
// temprature 

T=25 // degree celcius 

//length of plate, l=hn 
1=1 //m 

//width of plate ,w=0.5m 

w = 0.5 / /m 

//angle of incidence ,theta=0 degree 
theta=0 //degree 

/ / Calculation 

//for water at 25 degree celcius , momentum 
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diffusivity , 

MD = 8.63 *(10 ~-7) // m"2/s 

//local Reynold no. 
x = 0.5 / / m 
Re=x*v/MD 

//from Eq. 11.39, the boundary layer thickness is 

t=5*x/(Re "0.5) 


// Results 

printf (”i) Boundary layer thickness is%f m\n”,t) 
//local drag coefficient 

//CD=local drag force per unit area (F)/kinetic 
energy per unit volume (KE) 

//F = 0.332*rho *v"2* Re "0.5 and KE= 0.5*rho*v~2 

CD = 0.332 * v"2 *(Re"-0.5)/(0.5)*v"2 

printf (” Local drag coefficient is %f \n”,CD) 

//From eq 11.44, the drag force acting on one side 
of the plate is 
//kinetic viscocity 

mu=8.6*(10"-4) 

fd=0.664*mu*v*(l*v/MD)" 0.5*w 

//the total force acting on both sides of the plate 

tfd = 2* fd 

printf (”total drag force is %f N \n”,tfd) 

//shear stress at any point in the boundary layer 
//at a point in the boundary layer , 

x = 0.5 / / m 
y = t/2 

// n= blasius dimensionless variable 

n = y/(MD*x/v) "0.5 

//From table 11.1, at n = 2.5,f”(n) =0.218 
//shear stress= tau 
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fn = 0.218 / /f”(n) = fn 

tau = (mu*v*(v/(MD*x))~ 0.5)* fn 

printf (” Shear stress is %f N/n7 2”,tau) 


Scilab code Exa 11.2 air at 30 degree celcius 

//Example 11.2 
//Page no. 488 

//Calculate the thermal boundary layer thickness & 
//local heat transfer coefficient 0.75 m from the 
leading edge . 


//Variable 

Ts = 200 
Ta = 30 
V a = 8 
d=0.75 


declaration 

// C, temp . of air 

//C, temp .of surface 

//m/s, velocity of air 

//m, distant from leading edge 


II Calculation 

Tm=(Ts+Ta)/2 
mu=2.5*10"-5 
P = 0.69 
k=0.036 
Re=d*Va/mu 
t=5*d/(Re"0.5*P 


//C, Mean temp, of boundary layer 
//m" 2/s, viscosity 
//prndatl no. 

//W/m c, thermal conductivity 
/ / reynold no . 

(1/3)) //m, thermal 

boundary layer thickness 

printf (” Thermal boundary layer thickness is %f mm \n 
” ,t *10~3) 


N=(0.332*Re~(0.5)*P~(l/3)) //Nusslet no. 
h=k*N/d //heat 

transfer coefficent 

printf (”heat transfer coeff is %f W/m*2 C”,h) 
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Scilab code Exa 11.3 A thin metal plate 

//Example 11.3 
//Page No. 489 
// given 

//Free strean velocity (vl) and temp.(tl) on side 1 
vl=6 //m/s 

tl=150 //degree celcius 
//same on side 2 
v2 =3 //m/s 

t2=50 //degree celcius 
// distant 

x = 0.7 / / m 

//The plate temp, is assumed to be equal to the mean 
of the bulk air temp on the two sides of the 
plates 

T=100 //degree celcius 
//Side 1 

//mean air temp. 

tml=(T+t1)/2 

//From thermophysical properties: kinetic viscosity ( 
kv) ,Prandtl no.(P) , thermal conductivity (k) 

kvl=2.6*10 ~-5 //nP2/s 
P1=0.69 

kl=0.0336 //W/m degree celcius 
//Reynold no. 

Rel=x*vl/kvl 
//Nusslet no(Nl) 

a = l/3 

N1=0.332*(Re 1)"0.5*Pl~a 
hl=kl*Nl/x 

//Side 2 of the plate 

tm2 = (T + t2)/2 

/ / Similarly 
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kv2=2.076*(10)~-5 //nT2/s 
P2 =0.70 

k2=0.03 //W/m degree celcius 

Re2=x*v2/kv2 

N2 = 0.332*(Re2)“0.5*P2~a 

h2 = k2 *N2/x 

//overall heat transfer coeff. 

U=hl*h2/(hl+h2) 

//The local rate of heat exchange 

RH=U*(tl-t2) 

printf (” Local rate of heat exchange is %f W/m2\n\n” 

, RH) 

//the plate temp is given by 

TP=t2+(tl-t2)*U/h2 

printf (” Plate temperature is : %f Celsius \n”,TP) 


Scilab code Exa 11.4 calculate the temprature 


//Example 11.4 

//Calculate the temprature of the plate after 1 hour 
//if its initial temp , is 120 C 


//Given 

T1=120 

T2 =25 

Tm=(Tl+T2)/2 
rho =8880 

of plate 

//Properties of air at mean 

mu = 2.07*10 ~-5 

viscosity 

Pr=0.7 
k=0.03 

conductivity 


/ / C, initial temp 

//C, Final temp . 

//C, mean temp. 
//kg/m'3, density 

temp . 

//nP2/s , 

//Prandtl no . 
//W/m C, thermal 
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21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

1 

2 


1=0.4 

plate 

w = 0.3 
plate 
d=0.0254 
plate 
Vinf=1 

velocity 
Re=l*Vinf/mu 


//m, length of 
//m, width of 
/ /m, thickness of 
/ / m/ s , air 
// REynold no . 


// from eq . 11.90 (b) 

Nu=0.664*(Re)~(l/2)*(Pr)~(l/3) 
no . 

//Nu=1 *h/k 

h=Nu*k/l 

transfer coefficient 
//Rate of change of temp, is given 

A=2*l*w 
plate 
t=l*3600 
cp = 0.385*10 ~ 3 
specific heat 
m=l*w*d*rho 
plate 


//average Nusslet 

//W/hT 2 C, heat 
by 

//m "2. area of 

//s , time 
//j/kg K, 

//kg, mass of 


//—d/ dt (m* cp8dt )=A*hv* (Tl—T2) 

//appling the boundary condition 

T=(T1-T2)*exp(-A*h*t/(m*cp))+T2 

printf(”The temprature of plate after 1 hour is %f 
C” , round (T)) 


Scilab code Exa 11.5 Prandtl analogy 

//Example 11.5 
//Page no. 508 
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// given 

//Reynold no (Restriction factor ( f ), Prandlt no. (P) 

Re = 7.44*(10~4) 
f=0.00485 
P = 5.12 

x = P-l //assume 

//according to Von Karmen analogy 

N=((f/2)*Re*P)/(1+(5* sqrt (f/2))*(x+log(l+(5/6)*x))) 
printf (” Nusslet no is: %f \n”,N) 

// printf (” The prandtl analogy predicts Nu = 458.7”) 
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